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Introduction:  Prostate  cancer  (PCa)  is  the  most  common  cancer  in  men  in  several 
countries,  with  the  American  Cancer  Society  (ACS)  estimating  241,740  new  cases  of 
PCa  to  be  diagnosed  and  deaths  of  28,170  men  of  PCa  (1).  Due  to  its  prevalence  in  the 
male  population  as  well  as  its  unpredictable  clinical  course,  early  detection  and 
diagnosis  have  become  a  priority  for  many  health  care  professionals.  Another  method 
for  staging  prostate  cancer  is  through  imaging  techniques  including  ultrasound, 
computed  tomography  (CT),  and  magnetic  resonance  imaging  (MRI)  with  or  without  the 
help  of  dynamic  contrast  enhancement  modeling  (DCE-MRI),  diffusion  weighted  imaging 
(DWI),  and  magnetic  resonance  spectroscopy  (MRS)  (2-5).  MRS  is  a  powerful  tool  for 
exploring  the  cellular  chemistry  of  human  tissues  (3,5,6-11).  There  is  a  growing  body  of 
evidence  that  1H  MRS  may  contribute  to  the  clinical  evaluation  of  prostate  cancer  and 
also  for  evaluating  the  metabolic  alterations  due  to  therapy.  There  have  been  no  reports 
on  combining  two  spectral  dimensions  with  two-dimensional  (2D)  or  three  dimensional 
(3D)  spatial  encoding  applicable  to  prostate  cancer.  Acceleration  of  magnetic  resonance 
spectroscopic  imaging  (MRSI)  has  been  demonstrated  using  echo-planar  imaging 
techniques  (12-13).  Recently,  Schulte  et  al.  have  successfully  developed  an  algorithm 
called  prior-knowledge  fitting  (ProFit)  to  quantify  metabolite  concentrations  using  the 
JPRESS  spectra  recorded  using  a  Philips  3T  MRI  scanner  (14).  It  was  demonstrated 
that  metabolite  quantitation  of  JPRESS  spectra  with  ProFit  was  accurate,  robust  and 
yielding  generally  consistent  results,  both  in  vivo  and  in  vitro.  Their  results  suggest  that 
the  number  of  quantifiable  prostate  metabolites  can  be  increased  from  3-4  with  ID 
PRESS/LC-Model  to  more  than  10  with  JPRESS/ProFit  (15-16). 

Body: 

/)  Proposed  Task  1  (Months  1-6):  To  implement  a  multi-voxel  based  extension  of  the 
JPRESS  sequence,  in  which  two  spectral  encodings  will  be  combined  with  two  spatial 
encodings  using  the  new  Siemens  VB17a  platform.  This  four-dimensional  (4D)  data 
acquisition  scheme  will  be  accomplished  utilizing  the  EPI  approach  that  is  commonly 
used  for  spatial  encoding  in  MRI. 

Accomplished  during  September  29,  2011-October  28  2012:  The  4D  echo-planar  J- 
resolved  spectroscopic  imaging  (EP-JRESI)  sequence  as  shown  in  Fig.1  was 
successfully  compiled  using  the  Siemens  VB17a  compiler.  The  volume  of  interest  (VOI) 
was  localized  using  three  selective  radio-frequency  (RF)  pulses  similar  to  the  PRESS 
sequence  that  is  routinely  used  by  MR  researchers  around  the  world.  There  were  B0 
crusher  gradient  pulses  around  the  two  refocusing  180°  (tt)  RF  pulses.  The  EPI  read-out 
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enabled  readout  of  512  maximum  t2  spectral  points  and  16kx  spatial  points.  Remaining 
spectral  and  spatial  encodings  were  accomplished  using  incremented  At!  and  ky.  Total 
echo  time  (TE=  TE1+  TE2)  was  30  ms.  The  4D  raw  data  was  typically  512t2*64- 
100ti*16kx*16ky.  Since  the  4D  EP-JRESI  sequence  is  a  home-built  sequence,  there  is  no 
manufacturer  supplied  extraction  program  to  post-process  the  acquired  data.  Hence,  our 
group  had  to  develop  extraction  codes  using  the  MATLAB  library. 


TEi 


TE2 


Figure.1.  The  4D  EPJRESI  sequence  showing  two  spatial  encodings  (kx,ky)  and  two 
spectral  encodings  (t2,  ti).  Inherent  acceleration  of  kx  and  t2  is  enabled  here  by  the  echo- 
planar  imaging  (EPI)-based  read-out. 


ii)  Proposed  Task  2:  To  evaluate  the  EPI-based  JPRESS  using  a  prostate  phantom 
containing  several  metabolites  which  have  been  reported  in  prostate  tissues,  and  to 
optimize  the  EP-JJRESI  sequence  and  other  acquisition  parameters  using  the  phantom 

(Months  6-12). 


Accomplished  during  September  29,  2011-October  28  2012:  The  sequence  was 
tested  using  a  prostate  phantom  containing  10  different  metabolites  at  physiological 
concentrations  (pH  set  to  7.2).  A  500  ml  prostate  phantom  was  prepared  containing  the 
following  metabolites  at  physiological  concentrations  as  reported  in  healthy  human 
prostate  (15-16):  Citrate  (Cit,  50mM),  Creatine  (Cr,  5mM),  Choline  (Cho,1mM), 
Spermine  (Spm,  6mM),  myo-inositol  (ml,10mM),  Phosphocholine  (PCh,  2mM),  Taurine 
(Tau,  3mM),  Glutamate  (Glu,  4mM),  Glutamine  (Gin,  2.5mM)  and  Scyllo-lnositol  (si, 
0.8mM).  Shown  in  Fig.2A  is  an  axial  MRI  slice  image  showing  the  multi-voxel  grids  of 
MRSI  with  the  yellow  boundary  of  the  field  of  view  (FOV),  and  the  white  box  representing 
the  volume  of  interest  (VOI)  localized  by  the  PRESS  sequence  which  is  an  integral  part 
of  the  EP-JRESI  sequence.  The  following  parameters  were  used  for  acquiring  the  fully 
sampled  EP-JRESI  data:  TR/TE=1500/30ms,  16  phase  encodes  (ky),  32  read-out  points 
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(kx)  with  oversampling,  512  ±read-out  trains  resulting  in  512  pairs  of  complex  spectral 
points  in  the  2nd  spectral  dimension  (t2),  lOOh  increments  for  the  indirect  spectral 
dimension  and  one  average  per  encoding.  A  total  duration  of  40  minutes  was  necessary 
to  acquire  this  water-suppressed  EP-JRESI  data.  A  non-water-suppressed  EP-JRESI 
data  using  4  averages  with  only  one  ti  increment  was  used  for  eddy  current  and  phase 
correction  of  the  suppressed  data  (14).  After  apodization  and  Fourier  transformation  of 
this  4D  data,  the  reconstructed  multi-voxel  2D  J-resolved  spectra  (shown  in  C)  were 
overlaid  on  top  of  the  16x16  spatial  grids.  An  extracted  2D  J-resolved  spectrum  (3ml) 
around  the  center  of  the  VOI  is  shown  in  Fig.  2B. 

A  B 

Clt 

ml  Ch  Cr 


3.9  3.6  3.3  3.0  2.7  2.4  2.1  1.8  1.5 

F2  ( ppm ) 


C 


Figure  2.A.  An  MRI  axial  slice  image  of  the  prostate  phantom  showing  the  VOI 
localization  (white  box)  and  field  of  view  (FOV)  for  spatial  encoding  (yellow  grids).  B. 
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Extracted  3ml  2D  J-resolved  spectrum  from  the  center  of  the  phantom.  C.  Multi-voxel 
display  of  Cit  recordable  within  the  white  box  boundary;  the  peaks  outside  the  VOI  are 
due  to  Gibb's  ringing  or  bleed. 

As  mentioned  before,  Fig. 2  demonstrates  the  successful  implementation  of  the  4D  EP- 
JRESI  sequence  on  the  3T  MRI  scanner.  Significant  programming  efforts  were 
necessary  for  accomplishing  this  on  both  acquisition  and  post-processing  fronts. 

///>  Proposed  Task  3:  To  record  the  4D  EP-JRESI  spectra  in  the  peripheral,  central  and 
transition  zones  of  malignant,  benign  and  healthy  prostates.  (Months  6-12). 
Accomplished  during  September  29,  2011-October  28  2012:  A  32  yo  healthy  human 
subject  was  investigated  on  the  3T  MRI  scanner  using  the  quadrature  body  coil 
“transmit”  and  external  body  matrix  “receive”  coil  assembly.  The  following  parameters 
were  used  to  acquire  the  EP-JRESI  data:  TR/TE=1.5s/30ms,  2  averages,  512  t2, 
oversampled  32kx,  64  increments  along  the  indirect  spectral  (ti)  and  16  spatial  ky 
dimensions.  In  Fig.  3A,  an  axial  MRI  of  the  abdomen  is  shown  displaying  the  VOI 
covering  the  prostate  localized  by  the  PRESS  sequence  and  the  EP-JRESI  grids.  The 
multi-voxel  display  of  a  selected  region  is  shown  in  Fig.  3B.  The  2D  J-resolved  spectrum 
(2ml)  extracted  from  the  location  (*)  is  shown  in  Fig.  3C.  The  2D  diagonal  and  cross 
peaks  of  Cit  and  other  metabolites  are  visible.  As  shown  below,  the  endorectal  “receive” 
coil  is  expected  to  facilitate  at  least  one  order  of  magnitude  higher  sensitivity  than  that  of 
the  body  matrix  assembly  used  here. 


nd  Ch  ,57 


A  B  C 

Figure  3.  A.  An  axial  MRI  slice  of  the  32y.o.  healthy  male  subject  showing  the  VOI  and 
MRSI  grids.  B.  Expanded  multi-voxel  2D  spectra  showing  the  Cit  multiplets.  C.  The  2D 
J-resolved  spectrum  extracted  from  a  2ml  voxel  of  the  peripheral  zone. 


Using  the  endorectal  "receive"  coil,  another  4D  EP-JRESI  data  acquired  in  a  61 
year  old  PCa  patient  having  PSA  of  9.1  and  two  malignant  lesions  (GS3+4  in  the  right 
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base  and  GS4+3  in  the  right  mid  regions)  and  extracted  spectra  from  2  different 
locations  (1ml)  are  shown  in  Figure  4.  The  extracted  two  spectra  from  the  malignant 
masses  in  the  right  mid  (GS4+3)  (ii)  and  right  base  (iii)  regions,  and  iv)  healthy 
peripheral  zone. 


Ch 


F^O  (Hz) 


34  34  34  U  3.0  24  24  24  2.2  24 
F2  (ppm) 

\  iv)  Healthy 


Ch 


■  >• '  fan  &  o  (J3  F,=0  (Hz) 


Cit  j 

*  r 

T'f 


3.8  3.6  3.4  3.2  3.0  2.8  2.6  2.4  2.2  2.0 

F2  (pp*i) 


F,=0  (Hz) 


ii)  Cancer 


3.0  34  34  3.2  3.0  2.0  24  24  2.2  2.0 
F2  f.pP™) 

iii)  Cancer 


Figure  4.  i)  An  axial  T2W 
MRI  slice  showing  the 
multi-voxel  EP-JRESI 
grids;  extracted  2D 
JPRESS  spectra  of 
malignant  voxels  (ii)  and 
iii),  and  the  healthy  voxel 
in  the  peripheral  zone 
(iv). 


/V)  Proposed  Task  4;  To  develop,  evaluate  and  optimize  the  prior-knowledge  basis  set 
spectra  using  the  GAMMA-simulation  and  prostate  phantom  solutions  as  prior 
knowledge  for  the  multi-voxel  based  JPRESS  spectra  recorded  using  the  3T  MRI 
scanner  (Months  6-12). 

Accomplished  during  September  29, 2011-October  28  2012: 

Shown  in  Fig. 5  is  a  2D  J-Resolved  spectrum  recorded  in  a  27  y.o.  healthy  male. 
Presence  of  several  metabolites  is  clearly  marked  here.  More  than  10  basis-sets  were 
constructed  for  the  prior-knowledge  fitting  developed  using  the  GAMMA  library  (17). 
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Figure  6.  Prior-knowledge  J-resolved  spectra  for  the  prostate  metabolites:  Experimental 
(top),  fitted  using  ProFit  algorithm  (middle)  and  residual  (bottom). 


Shown  in  Fig. 6  are  2D  plots  showing  ProFit  quantitation  using  MATLAB  codes.  Using 
the  spectral  data  recorded  in  the  peripheral  zone  of  the  27y.o.  healthy  volunteer, 
following  ratios  were  quantified  (/Cr):  PCh:  0.09,  Cit:  7.36,  Spm:  8.55, 
Glycerlphosphocholine  (GPC):  0.483,  Cho:  0.01,  ml:3.79,  sl:0.012,  Tau:1.47,  Gin:  0.03, 
Glu:1 .52.  These  preliminary  results  show  the  feasibility  of  detecting  more  metabolites  in 
the  4D  EP-JRESI  spectrum  than  the  conventional  MRSI  showing  3-4  metabolites  only. 


Key  Research  Accomplishments 

•  Implementation  of  the  4D  EP-JRESI  sequence  on  the  Siemens  3T  and  1.5T  MRI 
scanners  after  the  C++  compilation  of  the  sequence  using  the  Siemens  IDEA 
VB17a  compiler. 

•  Prior-knowledge  basissets  have  been  developed  for  prostate  metabolites  such  as 
Cit,  choline  groups  containing  glycerylphosphocholine,  phosphocholine,  Cr,  Spm, 
ml,  Glu,  Gin,  si,  phosphoethanolamine  and  lactate  using  a  GAMMA  C++  library. 
Prostate  metabolite  quantitation  has  been  tested  using  the  ProFit  algorithm. 

•  Evaluation  of  the  4D  EP-JRESI  sequence  in  24  malignant  and  benign  prostate 
cancer  patients. 

•  Recent  progress  has  been  made  on  further  acceleration  of  the  4D  EP-JRESI 
sequence  using  non-uniform  undersampling  (NUS)  and  compressed  sensing 
(CS)  reconstruction.  This  will  enable  shortening  the  endorectal  spectral 
acquisition  and  reducing  the  patient  inconvenience  during  the  scan. 

•  Using  the  preliminary  results  obtained  using  this  IDEA  grant,  an  R01  application 
entitled  "EchoPlanar  J-resolved  Prostate  Cancer  Metabolite  Imaging  Using 
Compressed  Sensing"  was  submitted  to  the  National  Institute  of  health  (NIH)  in 
July201 2/October  2011.  In  the  October  2012  review,  my  revised  application 
received  a  non-fundable  32nd  percentile  score  only  as  the  National  Cancer 
Institute's  funding  is  currently  below  10th  percentile  only. 
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Reportable  Outcomes: 

A.  Peer-reviewed  Publications: 

1.  Nagarajan  R,  Margolis  D,  Raman  S,  et  al.  MR  Spectroscopic  Imaging  and  Diffusion- 
Weighted  Imaging  of  Prostate  Cancer  With  Gleason  Scores.  J  Magn  Reson  Imaging 
2012;36:697-703  [PMID:  22581787], 

2.  Nagarajan  R,  Margolis  D,  Raman  S,  et  al.  Correlation  of  Gleason  Scores  with 
Diffusion-Weighted  Imaging  Findings  of  Prostate  Cancer.  Advances  in  Urology  2012; 
2012:374805.  Epub  2011  Dec  15  [PMID:  22216026]. 

3.  McClure  TD,  Margolis  DJ,  Reiter  RE,  et  al.  Use  of  MR  imaging  to  determine 
preservation  of  the  neurovascular  bundles  at  robotic-assisted  laparoscopic 
prostatectomy.  Radiology.  2012  Mar;262(3):874-83.  Epub  2012  Jan  24  [PMID: 
22274837], 

4.  Furuyama  J,  Wilson  NE,  Burns  BL,  et  al.  Application  of  Compressed  Sensing  to 
Multidimensional  MR  Spectroscopic  Imaging  in  Human  Prostate.  Magn  Reson  Med 
2012;67:1499-1505  [PMID:  22505247],  “ 

B.  Presentations:  During  this  year,  the  first  abstract  entitled  “Novel  Multi-dimensional 
Magnetic  Resonance  Spectroscopic  Imaging:  Implementation  and  Pilot  Validation  in 
Prostate  and  Breast  Cancer  in  vivo”  was  submitted  and  presented  as  an  invited  talk  at 
the  17th  International  Biophysics  Congress  (IUPAB)  conference  in  Beijing,  China  (Oct. 30- 
Nov.3,  2011).  The  2nd  abstract  summarizing  the  implementation  of  the  4D  EP-JRESI 
sequence  and  evaluation  of  it  in  healthy  males  was  presented  at  the  2011  Radiological 
Society  of  Northern  America  (RSNA)  meeting  (Nov.-Dec.201 1)  in  Chicago,  IL.  2012. 
Third  abstract  entitled  "Accelerated  Multi-Voxel  Two-Dimensional  In  Vivo  Spectroscopy 
Using  Compressed  Sensing"  was  presented  at  the  53rd  Experimental  NMR  Conference 
(ENC)  meeting  in  Miami,  Florida  (April  15-20,  2012).  Fourth  Abstract  entitled 
"Accelerating  Echo-Planar  J-Resolved  Spectroscopy  of  the  Prostate  using  Compressed 
Sensing  in  a  Clinical  Setting"  was  presented  at  the  20th  International  Society  of 
Magnetic  Resonance  in  Medicine  (ISMRM)  meeting  in  Melbourne,  Australia  (May  5-11, 
2012). 

C.  Books:  None  on  Prostate  Cancer  Research  based. 

Conclusions:  First  few  months  of  the  1st  year  were  spent  in  getting  the  approval  from 
the  HSRRB  and  UCLA  IRB  offices.  The  4D  EP-JRESI  scanning  protocol  was 
successfully  implemented  on  the  3T  MRI  scanner.  After  optimizing  the  protocol  using 
phantom  solutions  containing  metabolites  and  corn  oil,  the  protocol  has  been 
successfully  tested  in  healthy  males,  and  malignant  and  BPH  patients.  We  will  continue 
to  recruit  20  malignant  and  5  benign  prostate  cancer  patients,  and  5  healthy  males 
during  the  next  year. 
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Purpose:  To  investigate  functional  changes  in  prostate 
cancer  patients  with  three  pathologically  proven  different 
Gleason  scores  (GS)  (3+3,  3+4,  and  4+3)  using  magnetic 
resonance  spectroscopic  imaging  (MRSI)  and  diffusion- 
weighted  imaging  (DWI). 

Materials  and  Methods:  In  this  study  MRSI  and  DWI 
data  were  acquired  in  4 1  prostate  cancer  patients  using  a 
1.5T  MRI  scanner  with  a  body  matrix  combined  with  an 
endorectal  coil.  The  metabolite  ratios  of  (Cho+Cr)/Cit 
were  calculated  from  the  peak  integrals  of  total  choline 
(Cho),  creatine  (Cr),  and  citrate  (Cit)  in  MRSI.  Apparent 
diffusion  coefficient  (ADC)  values  were  derived  from  DWI 
for  three  groups  of  Gleason  scores.  The  sensitivity  and 
specificity  of  MRSI  and  DWI  in  patients  were  calculated 
using  receiver  operating  characteristic  curve  (ROC) 
analysis. 

Results:  The  mean  and  standard  deviation  of  (Cho+Cr)/ 
Cit  ratios  of  GS  3+3,  GS  3+4,  and  GS  4+3  were:  0.44  ± 
0.02,  0.56  ±  0.06,  and  0.88  ±  0.11,  respectively.  For  the 
DWI,  the  mean  and  standard  deviation  of  ADC  values  in 
GS  3+3,  GS  3+4,  and  GS  4+3  were:  1.13  ±  0.11,  0.97  ± 
0.10,  and  0.83  ±  0.08  mm2 /sec,  respectively.  Statistical 
significances  were  observed  between  the  GS  and  metabo¬ 
lite  ratio  as  well  as  ADC  values  and  GS. 

Conclusion:  Combined  MRSI  and  DWI  helps  identify  the 
presence  and  the  proportion  of  aggressive  cancer  (ie, 
Gleason  grade  4)  that  might  not  be  apparent  on  biopsy 
sampling.  This  information  can  guide  subsequent 
rebiopsy  management,  especially  for  active  surveillance 
programs. 
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PROSTATE  CANCER  (PCa)  is  the  most  common  can¬ 
cer  in  men  and  is  the  second  leading  cause  of  cancer 
death  in  American  men,  behind  only  lung  cancer. 
Prostate  cancer  affects  over  200,000  men  each  year 
and  causes  27,000  deaths  annually  in  the  United 
States  (1).  Early  PCa  detection  is  the  key  to  successful 
cancer  treatment.  To  date,  the  suspicion  of  PCa  is 
mainly  based  on  elevated  serum  prostate-specific 
antigen  (PSA)  or  abnormal  digital  rectal  examination 
(DRE)  of  the  prostate,  both  with  well-known  limita¬ 
tions  (2,3).  Systematic  but  untargeted  transrectal 
ultrasound  (TRUS) -guided  biopsy  is  the  standard 
method  for  histological  diagnosis  of  PCa  (4) .  However, 
untargeted  systematic  biopsies  have  several  disadvan¬ 
tages:  they  miss  cancer  in  up  to  35%  of  cases  (4), 
multiple  foci  of  cancer  are  present  in  more  than  85% 
of  cases  of  PCa  (5),  and  complications  may  occur 
because  of  unnecessary  biopsies.  Therefore,  a  nonin- 
vasive  diagnostic  method  would  clearly  be  beneficial. 

The  pathologic  grade,  or  aggressiveness,  of  PCa  is 
given  by  the  Gleason  score  (6).  The  Gleason  grading 
system  consists  of  two  digits,  the  primary  and  second¬ 
ary  pattern  in  the  tumor.  Gleason  patterns  are  num¬ 
bered  in  increasing  order  of  cellular  disorder  and  loss 
of  normal  glandular  architecture  from  1  to  5,  with  1 
being  the  least  aggressive  and  5  the  most  aggressive. 
A  Gleason  score  (GS)  of  6  or  less  is  considered  low 
grade  (ie,  all  pattern  3  or  lower);  7  (ie,  3+4  or  4+3), 
intermediate  grade;  and  8  through  10  (some  compo¬ 
nent  of  pattern  4  or  5),  high  grade.  A  tumor  with  a 
low  GS  typically  grows  slowly  enough  that  it  may  not 
pose  a  significant  threat  to  the  patient  in  his  lifetime. 

Increased  choline  (Cho)  and  reduced  citrate  (Cit) 
have  been  demonstrated  by  3D  magnetic  resonance 
spectroscopic  imaging  (MRSI)  of  PCa.  Since  creatine 
(Cr)  is  not  resolvable  from  Cho  in  the  prostate  at  1.5T, 
combined  Cho+Cr  is  used.  The  decrease  in  Cit  in  PCa 
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is  due  to  both  changes  in  cellular  function  (eg, 
increased  metabolism  and  consumption  of  tricarb¬ 
oxylic  acid  cycle  metabolites)  and  in  the  organization 
of  the  tissue,  with  loss  of  its  characteristic  ductal 
morphology  (7,8).  The  elevation  of  the  Cho  peak  in 
PC  a  is  associated  with  changes  in  cell  membrane  syn¬ 
thesis  and  degradation  that  occur  with  the  evolution 
of  human  cancer  (9,10).  MR  spectroscopy  offers 
improved  sensitivity  and  specificity  for  PCa  detection 
(11-13).  There  are  several  reports  during  the  last  dec¬ 
ade  demonstrating  the  application  of  diffusion- 
weighted  imaging  (DWI)  in  PCa  and  these  findings 
look  promising  (14-17).  It  is  straightforward  to  imple¬ 
ment,  with  shorter  image  acquisition  times  and  the 
availability  of  quantitative  data  using  apparent  diffu¬ 
sion  coefficient  (ADC)  values. 

Magnetic  resonance  methods  have  shown  promising 
results  in  the  diagnosis  of  PCa  (18,19).  T2-weighted 
(T2W)  MRI  can  depict  prostate  zonal  anatomy  (20,21). 
MRSI  provides  information  about  relative  concentra¬ 
tions  of  cellular  metabolites  in  the  prostate  (22,23). 
DWI  characterizes  the  tissue  structure  at  the  micro¬ 
scopic  level  (24). 

The  purpose  of  the  study  was  to  assess  the  effec¬ 
tiveness  of  the  combined  3D  MRSI  with  DWI  investiga¬ 
tion  of  PCa  patients  with  a  progressively  increasing 
proportion  of  aggressive  cancer  (GS  3+3,  GS  3+4, 
and  GS  4+3). 

MATERIALS  AND  METHODS 
Patients 

Between  February  2007  and  June  2008,  41  patients 
(GS  3+3  [n  =  12],  GS  3+4  [n  -  20]  and  GS  4+3  [n  = 
9])  who  underwent  radical  prostatectomy  were 
selected  for  the  endorectal  MR  study.  At  least  6  weeks 
were  recommended  between  the  biopsy  and  MRI  due 
to  the  interference  of  postbiopsy  hemorrhage  with  the 
images.  The  Investigational  Review  Board  approved 
our  retrospective  study  and  waived  the  informed  con¬ 
sent  requirement.  The  ages  of  the  patients  ranged 
from  47  to  75  years,  or,  by  GS  category,  GS:  3  +  3 
(mean  ±  standard  deviation  [SD],  60.1  ±  6.7  years),  3 
+  4  (mean  ±  SD,  58.1  ±  4.2  years),  and  4  +  3  (mean 
±  SD,  60.3  ±  3.9  years).  The  mean  PSA  values  for 
the  patients  with  GS  3+3,  GS  3+4,  and  GS  4+3  were 
5.1  ng/ml,  6.7  ng/ml,  and  7.3  ng/ml,  respectively. 

MRI/DWI 

The  MRI/MRSI  data  were  acquired  on  a  1.5T  Avanto- 
Tim  MRI /MRS  scanner  with  high-performance 
gradients  (Siemens  Medical  Solutions,  Erlangen, 
Germany).  A  combined  body  matrix  phased-array  coil 
assembly  and  an  endorectal  coil  were  used  in  the 
“receive”  mode  while  a  quadrature  body  coil  was  used 
for  “transmit.”  After  acquiring  the  Tx -weighted  imaging 
(repetition  time  /  echo  time  [TR/TE]  =  540/11  msec), 
thin-section  high-spatial-resolution  transverse,  sagit¬ 
tal,  and  coronal  T2-weighted  fast  spin-echo  images  of 
the  prostate  and  seminal  vesicles  were  obtained  with 
the  following  parameters:  3800/101  msec  (effective); 


echo  train  length,  32;  section  thickness,  3  mm;  inter¬ 
section  gap,  0  mm;  field  of  view  (FOV),  14  cm;  matrix, 
205  x  256. 

Axial  DWI  images  were  recorded  using  the  single¬ 
shot  echo  planar  imaging  technique  using  the  follow¬ 
ing  imaging  parameters:  TR/TE  2000/83  msec, 
27  cm  FOV,  4-mm  slice  thickness,  0  mm  intersection 
gap,  3  averages.  Isotropic  DWIs  were  obtained  using 
diffusion  gradients  with  three  b-values  (0,  50,  and 
400  sec/mm2)  along  three  directions  of  motion-prob¬ 
ing  gradients.  Due  to  the  use  of  lower  field,  low  b-val¬ 
ues  were  chosen  for  the  DWI.  Higher  b-values  are 
associated  with  an  increased  background  noise.  Also, 
this  protocol  was  based  on  the  earlier  protocol  devel¬ 
oped  for  general  body  imaging.  The  orientation  and 
location  of  these  images  were  prescribed  identically  to 
the  transverse  T2-weighted  prostate  images. 

MRSI 

MRSI  was  performed  in  all  patients  with  water-  and 
fat-suppressed  acquisition.  The  3D  MRSI  parameters 
were  as  follows:  TR/TE  =  700/120  msec;  6  averages; 
spectral  width,  1300  Hz;  number  of  complex  points, 
512;  FOV,  80  x  80  x  80  mm3;  12  x  12  x  12  phase¬ 
encoding  steps  to  accomplish  spatial  encoding  along 
three  dimensions.  A  point- resolved  spectroscopy 
sequence  (PRESS)  was  used  to  acquire  the  proton  (XH) 
MR  spectra  from  a  volume  of  interest  of  55  x  40  x  40 
mm3.  The  MRSI  voxel  volume  was  0.3  cm3.  Outer  vol¬ 
ume  lipid  suppression  was  achieved  using  eight  3 -cm 
slab  saturation  pulses  around  the  volume  of  interest. 
Water  and  fat  resonances  were  suppressed  using  two 
12.6-msec  dual-frequency  selective  MEGA  pulses  (25) 
with  B0  crusher  gradients.  The  total  acquisition  time 
was  ^14  minutes  including  the  voxel  prescription, 
shimming,  and  prescan  optimization.  The  full 
width  at  half  maximum  (FWHM)  of  the  water  peak 
was  ^12  Hz. 

Data  Analysis:  MRSI  and  DWI 

The  MRSI  postprocessing  software  on  the  MR  scanner 
provided  by  the  manufacturer  was  used  to  evaluate 
the  (Cho+Cr)/Cit  ratios.  Postprocessing  of  the  time 
domain  data  extracted  from  MRSI  multiple  voxels  con¬ 
sisted  of  zero  filling  from  512  to  1024  data  points, 
multiplication  by  a  Hamming  filter,  Fourier  transfor¬ 
mation,  and  phase  and  baseline  correction.  We 
obtained  the  integral  values  by  fitting  Gaussian  line 
shape  functions  to  the  resulting  absorption  spectra. 
For  further  analysis,  ratios  of  integrals  (Cho+Cr)/Cit 
were  used.  The  MRSI  was  overlaid  onto  a  T2W  image, 
and  Digital  Imaging  and  Communications  in  Medicine 
(DICOM)  position  information  was  used  to  find  the 
closest  match.  In  patients  with  histologically  proven 
prostate  carcinoma,  the  voxels  covering  the  tumorous 
lesions  from  the  peripheral  zone  (PZ)  were  selected 
and  indicated  as  tumor  voxels.  Voxels  were  considered 
sufficient  for  analysis  if  one  of  the  metabolites  had  a 
signal-to-noise  ratio  (SNR)  of  at  least  4:1. 

In  the  DWI,  the  mean  ADCs  and  SD  of  peripheral 
zone  cancer  was  calculated  from  regions  of  interest 
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Table  1 

Statistically  Significant  Mean  ADC  Value,  MR  Spectroscopic  Ratio, 
and  Tumor  Volume  of  Three  Gleason  Scores 


Gleason 

scores 

(Cho+Cr)/ 

Cit  ratio  mean 
(range) 

Mean  ADC 
(range) 
(mm2/sec) 

Tumor 

volume 

(cc) 

3+3 

0.44  (0.34-0.49) 

1.13  (0.91-1.33) 

1 .54  (0. 5-3.5) 

3+4 

0.56  (0.46-0.73) 

0.97  (0.70-1.11) 

1 .96  (0.5-4. 6) 

4+3 

0.88  (0.75-1 .05) 

0.83  (0.74-1.01) 

2.77  (0.5-5) 

(ROIs)  defined  by  the  histological  findings  of  the  radi¬ 
cal  prostatectomy  specimens  and  compared  to  the  T2- 
weighted  imaging.  The  MRI  sections  and  histological 
slices  were  matched  based  on  the  level  sextant  loca¬ 
tions,  anterior /posterior  and  central/peripheral.  Each 
sextant  was  biopsied  twice  for  12  to  14  cores.  The  ROI 
of  the  same  size  and  location  from  each  method  were 
drawn  manually  on  the  ADC  maps  and  DWI  by  a  sin¬ 
gle  reader.  When  the  ROIs  were  drawn,  care  was 
taken  to  exclude  both  the  neurovascular  bundle  and 
the  urethra  to  reduce  any  error  in  ADC  calculations. 
If  a  tumor  was  located  in  several  imaging  slices  of 
ADC  maps,  ADC  values  were  measured  on  each  image 
of  the  ADC  maps  and  the  average  was  calculated.  The 
histology  was  reviewed  by  an  experienced  dedicated 
genitourinary  pathologist. 

Statistical  Analysis 

A  one-way  analysis  of  variance  (ANOVA)  was  per¬ 
formed  to  assess  the  statistical  differences  between 
MRSI  values  and  Gleason  scores  and  between  ADC 
values  and  Gleason  scores.  Spearman’s  rank  correla¬ 
tion  coefficients  were  used  to  characterize  the  associa¬ 
tion  of  the  ADC  values  and  MRSI  ratios  for  tumors 
with  the  Gleason  score  and  pathologically  determined 
tumor  volume. 
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Figure  1.  Box-and -whisker  plot  shows  the  MR  spectroscopic 
ratios  of  three  Gleason  scores  (GS  3+3,  GS  3+4,  and  GS 
4+3). 


Receiver  operating  characteristic  (ROC)  curve  analy¬ 
ses  based  on  logistic  regression  models  were  per¬ 
formed  in  order  to  identify  the  optimal  cutoff  value  for 
ADC  and  metabolite  ratios  for  prediction  purposes  of 
three  different  Gleason  scores.  The  area  under  the 
curve  (AUC),  interpreted  as  the  average  value  of  sensi¬ 
tivity  for  all  possible  values  of  specificity,  was  taken 
as  a  criterion  for  the  success  of  the  ROC  analysis.  An 
area  of  0.50  implies  that  the  variable  adds  no  infor¬ 
mation,  whereas  an  area  of  1  implies  perfect  accu¬ 
racy.  Sensitivity,  specificity,  positive  predictive  value 
(PPV),  and  negative  predictive  value  (NPV)  as  well  as 
accuracy  were  reported  for  the  optimal  thresholds.  P 
values  <0.05  were  considered  statistically  significant. 
Statistical  analyses  were  performed  using  SPSS  soft¬ 
ware  (Statistical  Package  for  the  Social  Sciences,  v. 
18.0,  Chicago,  IL). 

RESULTS 

Table  1  shows  the  mean  ADC  values  and  MRSI  ratios 
for  the  three  different  Gleason  score  categories  (GS 
3+3,  GS  3+4,  and  GS  4+3).  Figures  1  and  2  show  the 
“box  and  whisker”  plots  of  MRSI  ratios  and  ADC  val¬ 
ues  of  three  Gleason  scores  respectively.  Figure  3 
shows  the  multivoxel  display  of  a  59 -year -old  PCa 
patient  showing  elevation  of  Cho  and  Cr  in  the  right 
base  with  corresponding  ADC  map  on  the  right  side. 
In  the  MRSI,  the  mean  and  SD  of  (Cho+Cr)/Cit  ratio 
of  PCa  patients  with  GS  of  3+3  in  the  peripheral  zone 
was  0.44  ±  0.02.  For  GS  3+4  and  4+3  the  mean  and 
SD  of  (Cho+Cr)/Cit  ratios  were  0.56  ±  0.06  and  0.88 
±  0.11,  respectively.  Statistical  significance  (P  <  0.05) 
was  found  between  the  Gleason  score  categories  and 
the  metabolite  ratio  (Cho+Cr/Cit). 
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Figure  2.  Box-and-whisker  plot  shows  the  ADC  map 
three  Gleason  scores  (GS  3+3,  GS  3+4,  and  GS  4+3). 


of 


700 


Nagarajan  et  al. 


The  PCa  patients  with  higher  Gleason  scores  tended 
to  show  higher  metabolic  ratios  on  3D  MRSI.  In  the 
DWI  of  patients  with  GS  3+3,  the  mean  and  SD  of  the 
ADC  values  were  1.13  ±  0.11  mm2/sec  (35  ROIs).  For 
the  patients  with  GS  3+4  and  4+3,  the  mean  and  SD 
of  the  ADC  values  were  0.97  ±  0.10  mm2/sec  (48 
ROIs)  and  0.83  ±  0.08  mm2/sec  (27  ROIs),  respec¬ 
tively.  The  ROIs  were  drawn  to  encompass  the  largest 
area  of  PCa  without  including  the  tumor  margins,  the 
prostate  capsule,  or  the  urethra.  The  areas  of  ROIs 
ranged  from  4  to  10  mm2.  The  wide  range  of  areas 
was  due  to  the  wide  range  of  prostate  carcinomas. 
Measurements  were  conducted  through  a  circular  ROI 
on  lesions.  Statistical  significance  (P  <  0.05)  was 
observed  between  the  ADC  values  and  tumor  aggres¬ 
siveness.  The  ADC  intergroup  variability  between  the 
three  different  GS  falls  between  8%-10%  and  MRSI 
ratios  were  6%-13%.  The  intravariability  of  ADC  and 
MRSI  were  6%  and  9%,  respectively.  Although  there 
was  an  overlap  among  the  three  GS,  statistical  signifi¬ 
cances  were  observed  in  metabolites  ratio  and  ADC 
values.  The  reason  for  the  overlap  may  possibly  be 
due  to  heterogeneity  of  the  PCa. 

The  results  of  Spearman  correlation  analysis  are 
shown  in  Table  2.  We  found  a  strong  and  statistically 
significant  correlation  between  each  of  ADC  values 
and  MRSI  ratios  with  Gleason  score  (r  =  —0.665,  P  < 
0.001;  and  r=  0.902,  P  <  0.0001,  respectively).  There 
was  a  weak,  but  statistically  significant  correlation, 


Table  2 

Spearman  Rank  Correlation  to  Evaluate  Association  of  Each  of  the 
Measures  With  Gleason  Scores  and  Tumor  Volume 


Gleason  score 

Tumor  volume 

Measure 

Correlation 

P 

Correlation 

P 

ADC 

-0.66** 

<0.001 

0.36* 

.018 

(Cho+Cr)/Cit 

0.90** 

<0.001 

0.69** 

<0.001 

^Statistically  significant  at  the  0.001  level. 
^Statistically  significant  at  the  0.05  level. 


Figure  3.  T2-weighted  MRI 
with  MRSI  display  of  a  59- 
year-old  prostate  cancer 
patient  with  corresponding 
abnormal  ADC  (long  arrow) 
map. 

between  ADC  and  tumor  volume  (r  =  0.366,  P 
0.01 8)  and  a  strong  statistically  significant  correlation 
between  MRSI  ratio  and  tumor  volume  (r  —  0.694,  P  < 
0.001).  In  addition,  there  was  a  moderate  correlation 
that  reached  statistical  significance  between  tumor 
volumes  and  Gleason  scores  [r  =  0.360,  P  =  0.021). 

Figure  4  shows  the  comparison  of  ROC  curves  of 
ADC  and  MRSI  ratios  of  PCa  patients  with  Gleason 
scores  (Fig.  4A)  3+3  and  3+4,  (Fig.  4B)  3+3  and  4+3, 
(Fig.  4C)  3+4  and  4+3.  ROC  curve  analyses  for  differ¬ 
entiating  GS  3+3  and  3+4  suggest  an  optimal  cutoff 
value  of  1.05  x  10  3  mm2/s  for  ADC  value.  This 
implies  a  proportion  of  correctly  identified  true  GS 
3+4  cases  (sensitivity)  of  82.6%,  a  proportion  of  cor¬ 
rectly  identified  true  GS  3+3  cases  (specificity)  of 
68.2%,  a  proportion  of  true  GS  3+4  cases  (PPV)  of 
73.1%,  and  a  proportion  of  true  GS  3+3  cases  (NPV) 
of  78.9%.  The  AUC  was  84%  with  an  accuracy  of 
75.6%.  The  cutoff  values  for  the  MRSI  ratios 
(Cho+Cr)/Cit  were  0.49  with  sensitivity,  specificity 
and  accuracy  96.3%,  90.5%,  and  93.8,  respectively. 
The  AUC  of  the  MRSI  ratio  was  the  largest  (98.2%)  in 
the  differentiation  between  GS  3+3  from  3+4.  In  the 
evaluation  of  GS  3+3  versus  3+4  the  MRSI  ratios  had 
an  AUC  of  100%,  indicating  excellent  discrimination 
compared  to  AUC  of  98.6%  for  ADC.  Similarly,  AUC 
for  MRSI  ratios  for  distinguishing  GS  3+4  and  4+3 
was  100%  compared  to  AUC  of  87%  for  ADC.  Detailed 
results  of  the  ROC  curve  analyses  are  given  in  Table 
3,  giving  the  sensitivity,  specificity,  PPV,  NPV,  AUC, 
and  accuracy  for  classifying  the  cases  as  GS  3+3  ver¬ 
sus  3+4,  GS  3+3  versus  4+3,  and  as  GS  3+4  versus 
4+3,  respectively.  For  all  cases  the  addition  of  ADC 
with  MRSI  ratio  did  not  improve  the  classification  any 
further. 


DISCUSSION 

In  the  present  study  we  investigated  MRSI  and  DWI  of 
41  patients  with  pathologic  correlation  and  confirmed 
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—  ADC 

—  (Cho+Cr)/Cit 


Figure  4.  Graph  shows  ROC  curves  of  ADC  and  MRSI  ratios  of  prostate  cancer  patients  for  differentiation  of  GS  (a)  3+3 
from  3+4,  (b)  3+3  from  4+3,  and  (c)  3+4  from  4+3. 


that  the  ratios  of  (Cho+Cr)/Cit  in  the  lesion  are  posi¬ 
tively  correlated  with  the  lesion  Gleason  scores,  with 
an  increase  of  the  Cho  level  and  decrease  of  the  Cit 
level  with  rising  cancer  aggressiveness.  This  finding  is 
likely  clinically  relevant,  since  the  surgical  Gleason 
score  is  used  widely  in  models  that  are  used  to  predict 
outcome.  Combined  MRI  and  MRSI  aids  in  tumor 
localization  within  the  peripheral  zone  (26-28),  and 
preliminary  reports  of  MRSI  assessment  of  tumor 
aggressiveness  have  shown  promise  (29,30).  In  the 
last  decade,  many  publications  have  reported  the 
advantages  of  MRI  and  MRSI  combining  metabolic  in¬ 
formation  directly  on  the  corresponding  anatomic  dis¬ 
play  (31,32).  The  correlation  between  tumor  volume 
and  grade  agrees  with  published  reports  (33,34). 


A  significant  reduction  of  ADC  values  of  PC  a  is  well 
documented  (14-16,35-39)  and  is  further  supported 
by  our  investigation.  In  our  investigation  of  the  rela¬ 
tionship  between  the  diffusion  characteristics  and 
aggressiveness  of  PCa,  we  found  that  lower  ADC  val¬ 
ues  were  associated  with  higher  Gleason  scores  and 
that  ADC  values  may  help  to  differentiate  between 
low-risk  (GS  6)  and  intermediate-risk  (GS  7)  prostate 
cancer  if  the  tumor  is  visible  on  DWI.  DWI  as  a  func¬ 
tional  imaging  technique  can  measure  the  mobility  of 
water  within  tissues  in  addition  to  depicting  tumor 
size  and  shape.  DWI  may  detect  and  localize  PCa 
before  the  surgery  and  then  may  provide  qualitative 
or  quantitative  information  for  measuring  ADC  in 
patients  with  PCa  after  surgery.  We  found  that  mean 


Table  3 

Measures  of  Sensitivity,  Specificity,  PPV,  NPV,  and  Accuracy  of  ADC  and  MRSI  in  Discrimination  of  GS  3+3  From  3+4,  GS  3+3  From 
4+3,  and  GS  3+4  From  4+3  With  ROC  Curve  Analysis 


Parameter 

CV 

Sensitivity,  % 

Specificity,  % 

PPV,  % 

NPV,  % 

AUC,  % 

Accuracy,  % 

ADC 

GS  3+3  vs.  GS  3+4 

1.05 

82.6 

68.2 

73.1 

78.9 

84 

75.6 

GS  3+3  vs.  GS  4+3 

0.94 

92.3 

95.5 

92.3 

95.5 

98.6 

94.3 

GS  3+4  vs.  GS  4+3 

0.85 

69.2 

82.6 

69.2 

82.6 

87 

77.8 

(Cho+Cr)/Cit 

GS  3+3  vs.  GS  3+4 

0.49 

96.3 

90.5 

92.9 

95 

98.2 

93.8 

GS  3+3  vs.  GS  4+3 

0.75 

100 

100 

100 

100 

100 

100 

GS  3+4  vs.  GS  4+3 

0.75 

100 

100 

100 

100 

100 

100 

CV:  cutoff  value;  PPV:  positive  predictive  value;  NPV:  negative  predictive  value;  AUC:  area  under  the  curve. 
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ADC  values  were  significantly  lower  and  mean  metab¬ 
olite  ratios  values  were  significantly  higher  for  malig¬ 
nant  PZ  tissue  with  increasing  GS.  The  decrease  in 
ADC  values  in  malignant  tissue  is  attributed  to  histo¬ 
pathologic  characteristics,  including  hypercellularity, 
enlargement  of  nuclei,  hyperchromatisms,  and  angu¬ 
lation  of  the  nuclear  contour  (40),  which  result  in  a 
reduction  of  diffusional  displacement  of  water  mole¬ 
cules.  This  study  demonstrates  the  potentials  of  a 
combination  of  DWI,  MRSI,  and  tumor  volume  of 
patients  with  three  different  Gleason  scores  using 
ROC  curve  analysis. 

We  have  not  examined  the  correlation  between  bi¬ 
opsy  and  prostatectomy  Gleason  score  since  we  do 
not  have  the  biopsy  report  for  12  out  of  41  patients. 
In  13  out  of  29  (45%)  patients  the  biopsy  and  prosta¬ 
tectomy  Gleason  scores  were  the  same.  For  the 
remaining  16  patients,  five  patients  with  GS  3+3, 
seven  patients  with  GS  3+4,  and  four  patients  with 
GS  4+3  whose  biopsy  Gleason  scores  were  different 
from  prostatectomy.  The  distribution  of  metabolite 
values  for  GS  4+3  is  larger  than  the  other  group 
because  of  a  smaller  number  of  subjects.  In  this  study 
we  have  not  correlated  the  ADC  values  with  histology 
data.  We  observed  negative  correlation  between  Glea¬ 
son  scores  and  ADC  values,  and  positive  correlations 
between  Gleason  scores  and  MRSI  ratios,  and 
between  Gleason  scores  and  tumor  volumes. 

Most  of  the  studies  to  date  that  have  assessed  the 
utility  of  DWI  in  PCa  evaluated  lesions  within  the 
peripheral  zone  (PZ).  Several  studies  (41-44)  have 
suggested  improvement  in  tumor  detection  over  T2W, 
either  using  DWI  alone,  or  in  combination  with  T2W. 
Results  for  comparison  between  DWI  and  MRS  have 
been  mixed.  Chen  et  al  (45)  found  that  DWI  (82%)  and 
MRS  (84%)  had  similar  sensitivity  but  DWI  had  lower 
specificity  (82%)  compared  with  MRS  (98%).  In  con¬ 
trast,  Reinsberg  et  al  (46)  instead  showed  slightly 
lower  sensitivity  and  specificity  rates  of  DWI  com¬ 
pared  with  that  of  MRS,  particularly  when  the  per¬ 
centage  of  tumor  within  the  voxel  was  smaller. 

This  study  has  a  few  limitations.  In  the  DWI  and 
MRSI  in  PCa,  we  included  only  tumors  in  the  PZ, 
because  heterogeneity  in  the  hyperplastic  transition 
zone  of  older  men  limits  tumor  depiction.  However, 
the  majority  of  tumors  are  found  in  the  PZ.  Our  recent 
work  has  focused  on  further  investigation  using  the 
increased  sensitivity  advantages  of  the  3T  MRI  scan¬ 
ner  and  the  results  will  be  reported  later. 

In  conclusion,  combined  MRSI  and  DWI  can  help  in 
the  discrimination  of  intermediate  Gleason  grade 
tumors  from  low  Gleason  grade  tumors  with  histopa- 
thology  as  the  standard  of  reference  (“gold  standard”). 
This  information  can  guide  subsequent  definitive 
management  and  help  to  optimize  active  surveillance 
programs. 
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The  purpose  of  our  study  was  to  compare  the  apparent  diffusion  coefficient  (ADC)  derived  from  diffusion -weighted  imaging 
(DWI)  of  prostate  cancer  (PCa)  patients  with  three  classes  of  pathological  Gleason  scores  (GS).  Patients  whose  GS  met  these 
criteria  (GS  3  +  3,  GS  3  +  4,  and  GS  4  +  3)  were  included  in  this  study.  The  DWI  was  performed  using  b  values  of  0,  50,  and 
400  s/mm2  in  44  patients  using  an  endorectal  coil  on  a  1.5T  MRI  scanner.  The  apparent  diffusion  coefficient  (ADC)  values  were 
calculated  from  the  DWI  data  of  patients  with  three  different  Gleason  scores.  In  patients  with  a  high-grade  Gleason  score  (4  +  3), 
the  ADC  values  were  lower  in  the  peripheral  gland  tissue,  pathologically  determined  as  tumor  compared  to  low  grade  (3  +  3  and 
3  +  4).  The  mean  and  standard  deviation  of  the  ADC  values  for  patients  with  GS  3  +  3,  GS  3  +  4,  and  GS  4  +  3  were  1.135  ±  0. 119, 
0.976  ±  0.103  and  0.831  ±  0.087  mm2/sec.  The  ADC  values  were  statistically  significant  (P  <  0.05)  between  the  three  different 
scores  with  a  trend  of  decreasing  ADC  values  with  increasing  Gleason  scores  by  one-way  ANOVA  method.  This  study  shows  that 
the  DWI-derived  ADC  values  may  help  differentiate  aggressive  from  low-grade  PCa. 


1.  Introduction 

Prostate  cancer  (PCa)  is  the  most  common  malignancy 
among  men  in  the  USA,  with  an  estimated  217,730  new  cases 
and  32,050  PCa-related  deaths  in  2010  [1].  The  incidence 
of  PCa  increases  with  age,  and  it  is  very  uncommon  in 
men  younger  than  50  years  old.  With  greater  longevity  and 
increased  awareness  of  the  disease  leading  to  more  men 
requesting  screening,  it  is  to  be  expected  that  there  will  be  an 
increase  in  the  number  of  patients  diagnosed  with  PCa  in  the 
future.  Most  men  diagnosed  with  PCa  ultimately  survive  the 
disease  and  die  of  other  causes.  The  overall  5 -year  survival 
rate  is  99%  for  all  stages,  but  only  34%  when  there  are 
distant  metastases  [2].  The  aim  of  PCa  management  is  to 
identify,  treat,  and  cure  patients  with  aggressive  disease  that 
may  prove  fatal  but  to  avoid  overtreating  those  in  whom 
the  disease  is  unlikely  to  be  life  threatening.  Most  patients 
diagnosed  with  PCa  have  localized  disease  confined  to  the 


prostate.  A  small  number  with  high-grade  tumors  will 
progress  to  develop  local,  extracapsular  tumor  extension  and 
distant  metastases. 

Prostate  tumors  are  graded  according  to  their  patholog¬ 
ical  appearance  with  a  Gleason  score  (GS),  which  represents 
the  sum  of  the  dominant  and  subdominant  histological 
patterns  (grades).  High  GSs  indicate  aggressive  tumors  with 
increased  potential  for  local  and  distant  spread;  Gleason 
grading  has  been  shown  to  provide  a  spectrum  of  risk  for 
all  patients  [3].  Magnetic  resonance  imaging  (MRI)  provides 
incremental  value  to  clinical  findings  in  staging  patients 
of  intermediate  risk.  For  example,  organ-confined  disease 
implies  that  the  patient  may  benefit  from  local  therapy  such 
as  surgery  [4].  MRI  is  more  accurate  than  either  digital 
rectal  examination  (DRE)  or  transrectal  ultrasound  (TRUS) 
biopsy  in  preoperative  anatomical  localization  of  PCa  [5]. 
The  sensitivity  and  specificity  of  T2 -weighted  imaging  for 
PCa  vary  widely  due  to  differences  in  imaging  techniques, 
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reference  standards,  criteria  for  defining  disease  involvement 
on  MRI,  and  interobserver  variability  [6].  In  a  meta-analysis 
by  Sonnad  et  al.  T2  -weighted  imaging  showed  a  maximum 
joint  sensitivity  and  specificity  rate  of  74%  for  the  staging 
of  PCa  [7].  In  T2 -weighted  imaging,  regions  of  PCa  show 
decreased  signal  intensity  relative  to  normal  peripheral  (PZ) 
tissue  because  of  increased  cell  density  and  a  loss  of  prostatic 
ducts  [8].  This  finding  is  nonspecific,  however,  because  other 
diseases  such  as  prostatitis  or  hyperplasia  can  also  cause  low 
signal  intensity  in  T2-weighted  imaging  [9-12]. 

Diffusion-weighted  imaging  (DWI)  is  another  MR-based 
technique  that  probes  functional  characteristics  of  tissues. 
The  clinical  success  of  DWI  has  led  to  a  broadening  ap¬ 
plication  in  the  prostate  gland.  Rapid  changes  in  diffusion 
properties  can  be  identified  by  calculating  the  apparent 
diffusion  coefficient  (ADC).  Dickinson  et  al.  [13]  showed  the 
standardizing  multiparametric  magnetic  resonance  imaging 
(mpMRI)  for  PCa  detection,  localization,  and  characteri¬ 
zation.  The  use  of  DWI  as  a  tool  for  the  evaluation  and 
management  of  prostatic  cancer  has  grown  steadily  in  the 
past  decades  [  14-16] .  The  purpose  of  the  study  was  to  record 
DWI  and  to  compare  ADC  values  derived  from  DWI  in  PCa 
patients  with  three  different  Gleason  scores  (3  +  3,  3  +  4,  and 
4  +  3). 

2.  Materials  and  Methods 

A  total  of  44  clinically  localized  PCa  patients  who  underwent 
radical  retropubic  prostatectomy  between  January,  2007  and 
May,  2008  were  selected  for  this  study.  The  entire  protocol 
was  approved  by  the  institutional  review  board  (IRB),  and 
an  informed  consent  was  obtained  from  each  human  subject. 
The  ages  of  the  patients  ranged  from  47  to  75  years,  and  the 
patients  fell  into  three  different  groups  by  surgery  GS:  3  +  3 
(mean  ±  SD,  60.1  ±  6.7  years),  3  +  4  (mean  ±  SD,  58.1  ±  4.2 
years),  and  4  +  3  (mean  ±  SD,  60.3  ±3.9  years).  The  mean 
prostate- specific  antigen  (PSA)  value  for  the  patients  with  GS 
3  +  3,  GS  3  +  4,  and  GS  4  +  3,  respectively,  were  5.0  ng/mL, 
6.8  ng/mL,  and  7.4  ng/mL. 

All  patients  underwent  prostate  imaging  with  an  en¬ 
dorectal  coil  on  a  1.5  Tesla  Avanto-Tim  MRI  scanner  with 
high-performance  gradients  (Siemens  Medical  Solutions, 
Erlangen,  Germany).  Sequences  included  axial  turbo  spin- 
echo  (TSE)  T2 -weighted  imaging  through  the  prostate  and 
seminal  vesicles  (TR/TE,  3800/101  ms;  slice  thickness,  3  mm; 
no  interslice  gap;  field  of  view  (FOV),  140  mm,  matrix  205  X 
256,  slice  thickness  3  mm,  interslice  gap  0  mm,  echo-train 
number  32,  turbo  factor  13).  In  addition,  echo-planar 
diffusion-weighted  sequences  sensitized  in  three  orthogonal 
planes  (TR/TE  2000/83  ms,  bandwidth  1396  Hz  in  the  EPI 
frequency  direction)  with  b  values  of  0,  50,  400  s/mm2  were 
obtained  at  the  same  slice  positions  as  the  axial  T2 -weighted 
images.  Twelve  4-mm-thick  slices  with  no  interslice  gap 
(27  cm  FOV)  with  three  averages  provided  coverage  of  the 
prostate  with  an  image  acquisition  time  of  less  than  a  minute. 
Isotropic  ADC  maps  were  generated  with  the  system  software 
using  all  b  values  and  taking  an  average  value  for  the  two 
directions  of  diffusion  sensitization. 


MR  Images  were  initially  reviewed  without  clinical  infor¬ 
mation,  but  the  final  report  was  generated  after  the  clinical 
information  was  reviewed.  The  histology  was  reviewed  by 
an  experienced  pathologist.  The  edge  and  the  contour  char¬ 
acteristics  of  the  lesions  were  defined  using  the  same  slices 
on  which  regions -of- interest  (ROI)  analyses  were  performed. 
ROIs  were  drawn  independently  on  the  ADC  maps,  and 
differences  in  measurement  were  resolved  by  consensus. 
ROIs  were  drawn  in  the  tumor  PZ  of  all  the  44  PCa  patients. 
ADCs  were  calculated  for  all  slices  by 

a ln(Sl  —  So)  2/  a 

ADC  =  -  — - z——  mnr/s,  ( 1 ) 

(b i  -  b0) 

where  Si  is  the  signal  intensity  of  a  voxel  after  application  of 
a  diffusion  gradient  and  So  is  the  echo  magnitude  without 
diffusion  gradients  applied  (b  =  Os/mm2).  Diffusion 
sensitivity  is  determined  by  the  difference  between  b\  and  bo. 
If  multiple  tumors  were  present  in  the  peripheral  zone,  the 
average  ADC  value  was  recorded  for  each  lesion.  The  MRI 
sections  and  histological  slices  were  matched  on  the  basis  of 
the  sextant  level,  anterior/posterior,  and  peripheral/central 
(transitional). 

At  the  time  of  these  examinations,  other  sequences 
performed  as  part  of  the  routine  prostate  MRI  protocol  at 
our  institution  but  not  assessed  in  this  study  included  sagittal 
and  coronal  TSE  T2 -weighted  imaging  sequences  through  the 
prostate  and  seminal  vesicles. 

3.  Statistical  Analysis 

Statistical  analyses  were  performed  to  assess  the  statistical 
differences  between  ADC  values  for  the  three  different 
Gleason  scores  (GS  3  +  3,  GS  3  +  4,  and  GS  4  +  3)  using 
one-way  analysis  of  variance  (ANOVA)  with  SPSS  software 
package  assuming  parameters  were  normally  distributed.  A 
P  value  of  less  than  0.05  was  considered  to  indicate  a  statis¬ 
tically  significant  difference.  To  explore  for  any  relationship 
between  the  ADC  value,  tumor  volumes,  and  the  Gleason 
score,  Pearson  correlation  was  performed  on  the  data. 
Also,  analysis  of  covariance  (ANCOVA)  was  done  on  ADC 
values  of  different  Gleason  scores  with  tumor  volume  as  a 
covariate  to  see  its  effect. 

4.  Results 

The  patients  mean  and  standard  deviation  (SD)  of  age  and 
PSA  and  ADC  values  for  tumor  PZ  regions  of  three  Gleason 
scores  are  shown  in  Table  1.  Figure  1(a)  shows  the  T2- 
weighted  MRI  of  a  68-year-old  PCa  patient  with  GS  3  +  4 
and  Figure  1  (b),  corresponding  ADC  map  with  low  signal  on 
the  left  base  PZ.  Figure  2  illustrates  a  box  plot  of  ADC  values 
for  PCa  in  the  peripheral  zone  tissue  categorized  by  the  three 
Gleason  scores.  In  13  patients  with  GS  3  +  3,  the  (mean  ± 
SD)  ADC  value  was  1. 135  ±0.1 19  mm2/sec  using  32  ROIs.  In 
22  patients  with  GS  3  +  4,  the  (mean  ±  SD)  ADC  value  was 
0.976  ±  0.103  mm2/sec  using  52  ROIs.  In  9  patients  with  GS 
4  +  3,  the  (mean  ±  SD)  ADC  value  was  0.831  ±0.087  mm2/sec 
using  24  ROIs.  Although  a  statistically  significant  difference 
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(a)  (b) 


Figure  1:  (a)  T2-weighted  MRI  of  68  yo  prostate  cancer  patient  with  GS  3  +  4  and  (b)  corresponding  ADC  map  with  low  signal  on  the  left 
base  PZ. 


GS3  +  3  GS3  +  4  GS4  +  3 

Gleason  scores 

Figure  2:  Box- Whisker  plots  of  ADC  values  of  Gleason  score  (GS) 
3  +  3,  (GS)  3  +  4  and  (GS)  4  +  3  patients.  The  center  horizontal  line 
indicates  the  median. 


Table  1:  Patient  clinical  information  and  ADC  values. 


Gleason  scores 

Age  mean  ±  SD 

PSA  mean  ± 
SD  (ng/mL) 

ADC  values 
mean  ±  SD 
(mm2/sec) 

3  +  3  («=  13) 

60.1  ±6.7 

5.0  ±2.8 

1.135  ±0.119 

3  +  4  (m  =  22) 

58.1  ±4.2 

6.8  ±  1.7 

0.976  ±0.103 

4  +  3  («  =  9) 

60.3  ±  3.9 

7.4  ±  2.8 

0.831  ±  0.087 

P  value 

NS 

NS 

<0.05*** 

N  S — Nonsignificant 
***— Significant. 


existed  between  the  three  groups  (P  <  0.05),  a  certain  degree 
of  overlap  between  tissue  types  was  evident.  There  was  no 
statistical  significance  between  the  PSA  and  patients  ages 
with  three  different  Gleason  scores. 


We  did  not  have  the  biopsy  report  for  14  patients  out 
of  44  patients.  Out  of  which  15  patients  whose  biopsy 
and  prostatectomy  Gleason  scores  were  the  same.  For  the 
remaining  15  patients,  We  had  four  patients  with  GS  3  + 
3,  eight  patients  with  GS  3  +  4,  and  3  patients  with  GS  4  + 
3  whose  biopsy  Gleason  scores  were  different  from  prosta¬ 
tectomy.  FFence,  we  have  not  done  the  correlation  between 
biopsy  and  prostatectomy  Gleason  scores.  Out  of  44  patients, 
in  35  patients  (GS  3  +  3  [n  =  13),  GS  3  +  4  {n  =  13), 
and  GS  4  +  3  ( n  =  9))  tumors  were  detected  by  the  DWI 
technique.  Nine  patients  (GS  3  +  4  [n  =  9))  were  missed 
by  the  MRI  technique.  To  evaluate  the  association  between 
ADC  value,  tumor  volumes,  and  the  Gleason  grade,  all  the 
ADC  values  and  tumor  volumes  were  first  summarized  at 
the  individual  patient  level,  followed  by  applying  Pearson  s 
correlation  coefficients.  We  observed  negative  correlation 
between  Gleason  score  and  ADC  values  and  positive  cor¬ 
relation  between  Gleason  score  and  tumor  volume.  In  the 
ANCOVA  analysis,  the  results  were  statistically  significant 
(P  =  0.0001)  between  the  Gleason  score  and  ADC  values. 

5.  Discussion 

To  increase  the  accuracy  of  MRI,  a  number  of  authors  have 
used  special  techniques  to  study  a  particular  characteristic  of 
the  prostate  tumor  and  surrounding  tissues  such  as  dynamic 
contrast-enhanced  (DCE)  MRI  [17-19]  and  MR  spec¬ 
troscopy  (MRS)  [20-24].  MR  spectroscopy  requires  a  sub¬ 
stantially  longer  examination  time  than  DWI,  and,  addition¬ 
ally,  shimming  process  and  placement  of  saturation  bands 
during  the  examination  are  time  consuming.  For  evaluation 
of  MRS,  baseline  correction  and  phase  correction  have  to 
be  performed  in  some  cases. 

DWI  is  the  only  functional  imaging  technique  that  is 
able  to  assess  molecular  diffusion  in  vivo  and  provides  infor¬ 
mation  about  biophysical  properties  of  tissues  such  as  cell 
organization,  density,  and  microstructure  [25].  DWI  maybe 
helpful  in  differentiating  high-risk  patients  from  those  at  low 
and  intermediate  risks,  since  there  is  a  significant  correlation 
between  the  ADC  values  from  patients  with  three  different 
Gleason  scores.  The  patients  with  the  Gleason  score  of  4  +  3 
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have  a  higher  likelihood  of  biochemical  recurrence,  partic¬ 
ularly  for  the  increasing  proportion  of  patients  with  organ- 
confined  disease  after  radical  prostatectomy  than  those  with 
3  +  4  as  reported  by  Sakr  et  al.  [26] .  Also,  our  results  showed 
decreased  ADC  values  in  patients  with  GS  4  +  3  than  those 
with  GS  3  +  4  significantly  (P  <  0.05).  This  may  be  useful  to 
assess  the  aggressiveness  of  the  PCa.  The  sensitivity  of  DWI 
is  better  in  the  PZ  than  the  central  gland  [27].  DWI  has 
also  been  shown  to  be  helpful  in  the  identification  of  PCa 
in  patients  with  previous  negative  biopsies  and  persistently 
elevated  PSA  [28]. 

McNeal  et  al.  [29]  reported  that  65%  of  prostate  cancers 
arise  in  the  PZ  and  up  to  30%  arise  in  the  transition 
zone  (TZ).  The  transrectal  MRI  is  generally  considered  less 
specific  for  use  in  the  evaluation  of  TZ  cancers  because  of 
the  heterogeneously  low  T2  signal  intensity  in  normal  TZ 
[30],  and  the  presence  of  benign  prostatic  hyperplasia  in  TZ 
[31,  32].  The  endorectal  coil  offers  poor  signal  sensitivity 
when  it  comes  to  the  TZ.  Hosseinzadeh  and  Schwarz  [33] 
successfully  investigated  T2  relaxation  rates  and  diffusion- 
weighted  images  of  the  human  prostate  using  an  endorectal 
coil.  DWI  can  provide  valuable  cellular  information  about 
tissue  in  addition  to  the  conventional  Ti  and  T2 -weighted 
imaging  [34-38].  In  this  present  study,  ADC  values  show 
a  decreasing  trend  with  increasing  Gleason  scores.  The 
calculated  ADCs  for  cancer  in  the  PZ  were  consistent 
with  those  previously  reported  studies  [25,  39-41].  These 
findings  suggest  that  measurement  of  ADC  may  provide  an 
additional  feature  that  could  further  increase  the  specificity 
of  diagnosis  for  PCa.  The  variation  in  reported  ADC  values 
could  be  due  to  a  number  of  physiologic  factors  (e.g.,  age 
tumor  size)  as  well  as  technical  factors  (e.g.,  variations 
in  acquisition  parameters,  inhomogeneous  signal  reception 
using  the  endorectal  MRI  coil,  and  postprocessing  methods). 

Our  study  suffered  some  limitations.  We  analyzed  cancer 
only  in  the  PZ,  where  most  cancers  occur.  DWI  itself  also 
has  some  limitations  as  this  sequence  is  affected  by  magnetic 
susceptibility,  resulting  in  spatial  distortion  and  signal  loss. 
Moreover,  there  is  no  consensus  on  the  optimal  b  value  for 
DWI  of  the  prostate. 

In  conclusion,  this  study  shows  that  the  DWI  correlates 
with  pathological  Gleason  scores.  DWI -acquired  ADC  values 
are  a  very  potential  measure  to  delineate  prostate  carcinoma 
from  the  PZ  and  are  able  to  predict  the  presence  of  low, 
and  high-grade  components  in  PCa  with  great  accuracy.  The 
ADC  values  derived  from  the  1.5  T  diffusion-weighted  MRI 
demonstrate  tumor  aggressiveness  and  could  be  of  future 
use  in  treatment  decisions  and  in  patient  followup  in  active 
surveillance. 
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Use  of  MR  Imaging  to 
Determine  Preservation  of 
the  Neurovascular  Bundles  at 
Robotic-assisted  Laparoscopic 
Prostatectomy1 


Purpose:  To  determine  whether  findings  at  preoperative  endorec¬ 

tal  coil  magnetic  resonance  (MR)  imaging  influence  the 
decision  to  preserve  neurovascular  bundles  and  the  extent 
of  surgical  margins  in  robotic-assisted  laparoscopic  pros¬ 
tatectomy  (RALP). 

Materials  and  This  study  was  approved  by  the  investigational  review 
Methods:  board  and  was  compliant  with  the  HIPAA;  the  require¬ 

ment  to  obtain  informed  consent  was  waived.  The  au¬ 
thors  prospectively  evaluated  104  consecutive  men  with 
biopsy-proved  prostate  cancer  who  underwent  preop¬ 
erative  endorectal  coil  MR  imaging  of  the  prostate  and 
subsequent  RALP.  MR  imaging  was  performed  at  1.5  T 
between  January  2004  and  April  2008  and  included  T2- 
weighted  imaging  ( n  =  104) ,  diffusion-weighted  imaging  ( n 
=  88),  dynamic  contrast-enhanced  imaging  ( n  =  51),  and 
MR  spectroscopy  (n  =  91).  One  surgeon  determined  the 
planned  preoperative  extent  of  resection  bilaterally  on  the 
basis  of  clinical  information  and  then  again  after  review  of 
the  final  MR  imaging  report.  The  differences  in  the  surgi¬ 
cal  plan  before  and  after  review  of  the  MR  imaging  report 
were  determined  and  compared  with  the  actual  surgical 
and  pathologic  results  by  using  logistic  regression  analysis. 
Continuous  and  ranked  variables  underwent  Pearson  and 
Spearman  analysis. 


Results:  After  review  of  MR  imaging  results,  the  initial  surgical 

plan  was  changed  in  28  of  the  104  patients  (27%);  the  sur¬ 
gical  plan  was  changed  to  a  nerve-sparing  technique  in  17 
of  the  28  patients  (61%)  and  to  a  non-nerve-sparing  tech¬ 
nique  in  11  (39%).  Seven  of  the  104  patients  (6.7%)  had 
positive  surgical  margins.  In  patients  whose  surgical  plan 
was  changed  to  a  nerve-sparing  technique,  there  were  no 
positive  margins  on  the  side  of  the  prostate  with  a  change 
in  treatment  plan. 


Conclusion:  Preoperative  prostate  MR  imaging  data  changed  the 

decision  to  use  a  nerve-sparing  technique  during  RALP  in 
27%  of  patients  in  this  series. 
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McClure  et  al 


Prostate  cancer  is  the  most  com¬ 
mon  noncutaneous  cancer  di¬ 
agnosed  in  American  men  (1). 
One  standard  of  care  for  therapy  of 
organ-confined  cancers  has  been  the 
radical  retropubic  prostatectomy, 
which  carries  a  substantial  risk  of 
morbidity — including  incontinence 

and  impotence.  A  variety  of  newer 
“nerve-sparing”  open  and  laparoscopic 
techniques  have  been  developed,  pur¬ 
portedly  to  address  this  issue  (2),  in¬ 
cluding  robotic-assisted  laparoscopic 
prostatectomy  (RALP)  (3).  RALP 
provides  improved  visualization  of 
the  surgical  field  and  improved  instru¬ 
ment  control  when  compared  with 
other  techniques.  However,  surgeons 
performing  RALP  lack  tactile  (haptic) 
feedback,  upon  which  they  tradition¬ 
ally  rely  to  determine  the  extent  of  re¬ 
section  (4).  Lack  of  tactile  sensation 
may  compromise  a  surgeon’s  ability 
to  appropriately  evaluate  potential  in¬ 
volvement  of  the  neurovascular  bun¬ 
dles  by  prostatic  adenocarcinoma. 
Surgeons  must  balance  the  extent  of 
resection  proportionate  to  the  volume 
of  cancer  with  the  probability  of  ex- 
tracapsular  infiltration.  Overly  aggres¬ 
sive  resections  risk  impotence  and 
incontinence,  but  overly  conservative 
resections  risk  positive  surgical  mar¬ 
gins.  No  conventional  preoperative 
urological  technique  provides  infor¬ 
mation  to  supplant  haptic  feedback. 

Magnetic  resonance  (MR)  imag¬ 
ing  of  the  prostate  with  use  of  an  en¬ 
dorectal  coil  spatially  localizes  pros¬ 
tate  cancers,  especially  those  that  are 
of  a  higher  volume  and  higher  grade 


Advances  in  Knowledge 


■  MR  imaging  information  influ¬ 
enced  the  extent  of  surgery 
(nerve-sparing  vs  non-nerve¬ 
sparing  technique)  in  28  of  the 
104  patients  (27%)  who  under¬ 
went  robotic-assisted  laparo¬ 
scopic  prostatectomy  (RALP) . 

■  Adding  MR  imaging  data  to  the 
clinical  information  to  determine 
the  extent  of  surgery  did  not 
result  in  additional  positive  mar¬ 
gins  at  surgery  (0  of  28  patients) . 


(5-12).  Preoperative  knowledge  of 
tumor  location  and  pathologic  stage, 
that  is,  the  presence  of  extracapsular 
extension  or  seminal  vesicle  invasion, 
provides  the  robotic  surgeon  addi¬ 
tional  information  that  can  be  used 
to  plan  the  surgical  technique.  The 
purpose  of  this  study  was  to  deter¬ 
mine  whether  findings  at  preoperative 
endorectal  coil  MR  imaging  influence 
the  decision  to  preserve  neurovascu¬ 
lar  bundles  and  the  extent  of  surgical 
margins  in  RALP. 


Materials  and  Methods 


Patient  Selection 

The  study  was  approved  by  our  insti¬ 
tutional  review  board  and  was  com¬ 
pliant  with  the  Health  Insurance  Por¬ 
tability  and  Accountability  Act.  The 
institutional  review  board  determined 
that  the  requirement  to  obtain  spe¬ 
cific  consent  for  the  study  was  not 
necessary  because  treatment  was 
within  the  standard  of  care  (written 
informed  consent  was  obtained  ex¬ 
plicitly  for  MR  imaging  and  surgery, 
with  the  understanding  that  surgical 
extent  would  be  based  on  the  imaging 


Implications  for  Patient  Care 


■  Prostate  MR  imaging  data  com¬ 
plement  clinical  information  to 
help  surgeons  determine  whether 
to  perform  nerve-sparing  or 
non-nerve -sparing  RALP  in  27% 
of  patients. 

■  Some  patients  thought  to  require 
non-nerve -sparing  surgery  might 
safely  undergo  nerve-sparing  sur¬ 
gery  on  the  basis  of  MR  imaging 
findings  because  this  change  in 
surgical  plan  does  not  appear  to 
be  associated  with  an  increase  in 
positive  surgical  margins. 

■  MR  imaging  can  help  localize  the 
side  of  stage  T3  disease  so  that 
only  an  ipsilateral  neurovascular 
bundle  may  be  resected;  the  con¬ 
tralateral  side  may  safely  un¬ 
dergo  nerve-sparing  resection  to 
minimize  postoperative  inconti¬ 
nence  or  impotency. 


findings).  Retween  January  2004  and 
April  2008,  we  prospectively  identi¬ 
fied  105  consecutive  patients  who  met 
the  inclusion  criterion  of  MR  imaging 
of  the  prostate  with  an  endorectal  coil 
before  RALP  (Table  1);  one  patient 
who  underwent  open  prostatectomy 
rather  than  RALP  was  excluded.  Ri- 
opsy-proved  adenocarcinoma  of  the 
prostate  gland  was  diagnosed  in  all 
patients  before  MR  imaging. 

Demographics 

Patient  cancer  grade  and  stage  are 
summarized  in  Table  2.  Not  all  patients 
underwent  all  components  of  the  MR 
imaging  prostate  examination.  All  pa¬ 
tients  underwent  T2-weighted  MR  im¬ 
aging.  In  addition  to  T2-weighted  MR 
imaging,  16  patients  underwent  MR 
spectroscopy,  36  patients  underwent 
both  MR  spectroscopy  and  diffusion- 
weighted  imaging,  12  patients  under¬ 
went  both  diffusion-weighted  imaging 
and  dynamic  contrast  material-en¬ 
hanced  imaging,  and  one  patient  un¬ 
derwent  diffusion-weighted  imaging. 
Thirty-nine  patients  underwent  all 
four  components.  The  variability  of 
MR  components  occurred  for  two  rea¬ 
sons.  First,  spectroscopic  imaging  was 
not  covered  by  all  insurers;  therefore, 
not  all  patients  elected  to  pay  out 
of  pocket  for  spectroscopy.  Second, 
both  diffusion-weighted  and  dynamic 
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Table  1 


Patient  Demographics 

Characteristic 

Mean 

Median 

Maximum 

Minimum 

Standard  Deviation 

Age  (y) 

60.1 

59.7 

85 

41 

6.8 

Serum  PSA  level  (ng/dL) 

6.5 

5.3 

43.10 

1.01 

4.7 

Time  between  MR  imaging  and 
surgery  (wk) 

7.7 

4.9 

56 

0 

9.1 

Time  between  biopsy  and  MR 
imaging  (wk) 

12.7 

11.3 

46 

2 

6.9 

contrast-enhanced  imaging  became 
available  at  different  times  during  the 
course  of  the  study. 

MR  Imaging 

All  104  patients  underwent  imaging 
with  a  1 .5-T  unit  with  high-performance 
gradients  (Avanto;  Siemens  Medical  So¬ 
lutions,  Erlangen,  Germany)  by  using 
an  inflatable  endorectal  coil  (Medrad, 
Indianola,  Pa)  insufflated  with  60-90 
mL  of  room  air.  All  patients  under¬ 
went  sagittal,  coronal,  and  axial  oblique 
turbo  spin-echo  T2- weighted  imaging. 
In  88  of  the  104  patients  (85%),  echo- 
planar  diffusion- weighted  imaging  was 
performed  with  calculation  of  apparent 
diffusion  coefficients.  In  51  of  the  104 
patients  (49%),  dynamic  contrast-en¬ 
hanced  imaging  was  performed  by  using 
time-resolved  imaging  with  stochastic 
trajectories,  a  k-space-sharing  three- 
dimensional  gradient-recalled  echo 
sequence  (with  0.1  mg/kg  gadopen- 
tetate  dimeglumine  injected  at  2  mL/ 
sec  with  the  third  of  14  acquisitions; 
reconstructed  voxel  size  was  0.9  X  1.2 
X  1.5  mm).  MR  spectroscopy  was  per¬ 
formed  in  91  of  the  104  patients  (88%). 
An  echo  time  of  120  msec  was  used  for 
optimal  citrate  resonance.  Data  were 
processed  online  by  using  the  Siemens 
postprocessing  package.  Imaging  pa¬ 
rameters  are  given  in  Table  3. 

Image  Analysis 

All  images  were  analyzed  prospectively 
by  consensus  of  two  fellowship -trained 
genitourinary  radiologists  (S.S.R.  and 
D.J.A.M.,  with  12  and  5  years  of  expe¬ 
rience,  respectively),  who  were  initially 
blinded  to  clinical  data.  Criteria  for 
the  presence  of  cancer  on  T2-weighted 


images  included  round,  ovoid,  or  irreg¬ 
ular  areas  of  low  signal  intensity  on  T2- 
weighted  images  without  corresponding 
high  T1  signal  intensity,  which  might 
indicate  hemorrhage.  Criteria  for  ex- 
tracapsular  extension  at  MR  imaging 
that  have  been  found  to  be  specific  in 
previous  studies  included  the  presence 
of  low  signal  intensity  in  the  peripheral 
zone  of  the  prostate  with  irregular  bulg¬ 
ing  or  bowing  of  the  prostate  capsule 
at  T2-weighted  imaging,  disruption  of 
the  low- signal-intensity  periprostatic 
band  at  T2 -weighted  imaging,  or  di¬ 
rect  involvement  of  the  neurovascular 
bundle  and/or  obliteration  of  the  ret- 
roprostatic  angle  (5,11-13).  The  likeli¬ 
hood  of  extracapsular  extension,  as  a 
composite  score  based  on  T2-  and  dif¬ 
fusion-weighted  imaging,  dynamic  con¬ 
trast-enhanced  imaging,  and  MR  spec¬ 
troscopy  as  available,  was  prospectively 
reported  as  low,  intermediate,  or  high 
for  each  side,  with  intermediate  and 
high  cases  combined  for  binary  analysis 
because  this  was  the  consideration  of 
the  surgeon.  When  apparent  diffusion 
coefficient  maps,  dynamic  contrast-en¬ 
hanced  images,  and  MR  spectroscopic 
images  were  available,  these  were  also 
used  to  inform  the  degree  of  suspicion. 
Lesions  were  considered  suspicious  if 
at  least  one  of  the  following  findings 
was  present:  restricted  diffusion,  early 
or  intense  enhancement,  washout,  or 
elevated  choline -to -citrate  ratio.  All 
MR  spectroscopic  data  were  obtained 
and  analyzed  under  the  supervision  of 
a  physicist  with  more  than  15  years  of 
experience  in  endorectal  MR  imaging 
(M.A.T.).  Voxels  were  considered  suffi¬ 
cient  for  analysis  if  one  metabolite  had 
a  signal-to-noise  ratio  of  at  least  4:1. 


Voxels  were  considered  positive  for  tu¬ 
mor  when  the  ratio  of  choline  and  cre¬ 
atine  to  citrate  was  at  least  1:3  on  at 
least  four  contiguous  voxels. 

Imaging  results  were  evaluated,  and 
disease  stage  was  assigned  according  to 
the  American  Joint  Committee  on  Can¬ 
cer  Staging  of  Prostate  Cancer  (TNM 
system);  thus,  neurovascular  bundle  in¬ 
volvement  was  not  separately  recorded. 

Robotic  Prostatectomy  Algorithm 

RALP  was  performed  by  one  surgeon 
(R.E.R.,  with  more  than  6  years  of  ex¬ 
perience  with  this  technique)  with  use 
of  a  Da  Vinci  robot  (Intuitive  Surgical, 
Sunnyvale,  Calif).  The  surgeon  initially 
formulated  a  resection  plan  (bilateral 
nerve-sparing  surgery,  bilateral  non¬ 
nerve-sparing  surgery,  or  a  combina¬ 
tion  of  unilateral  non-nerve-sparing 
and  nerve-sparing  surgery)  blinded  to 
MR  results  on  the  basis  of  a  compos¬ 
ite  of  the  following  clinical  variables: 
results  of  digital  rectal  examination, 
results  of  sextant  prostate  biopsy 
(including  Gleason  score,  number  of 
positive  core  specimens,  percentage 
of  positive  core  specimens,  and  tumor 
length),  and  serum  prostate-specific 
antigen  (PSA)  level.  The  resection 
plan  was  reevaluated  after  review  of 
the  MR  imaging  results.  The  nerve¬ 
sparing  technique  was  reserved  for 
patients  without  intermediate  or  high 
probability  of  extracapsular  extension 
or  neurovascular  bundle  involvement. 
The  surgical  technique  was  document¬ 
ed  in  the  surgery  report. 

Surgical-Pathologic  Examination 

The  prostate  was  serially  sectioned  (not 
whole  mount)  from  base  to  apex  into 
different  levels  (depending  on  the  size 
of  the  prostate)  for  histologic  analysis 
and  labeled  as  right  or  left  and  anterior 
or  posterior  apex,  midgland,  and  base. 
The  right  and  left  seminal  vesicle  and 
vas  deferens  were  also  analyzed  sepa¬ 
rately.  Slides  were  reviewed  by  experi¬ 
enced  genitourinary  pathologists  who 
were  not  aware  of  imaging  findings  per 
standard  clinical  protocol.  All  reports 
were  reviewed  to  determine  the  pres¬ 
ence  of  extracapsular  extension  or  sem¬ 
inal  vesicle  invasion  and  to  compare 
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Table  2 


Biopsy  and  Surgical  Specimen  Gleason  Score  with  Surgery  Type 


Biopsy  Gleason  Score 

No.  of  Patients 

Surgical  Specimen  Gleason  Score 

Nerve-Sparing  Surgery 

3  +  3  =  6 

3  +  4  =  7 

4  +  3  =  7 

4  +  4  =  8 

Initial  Plan 

Actual  Surgery 

3  +  3  =  6 

60  (58) 

27 

29 

3 

1 

47 

48 

3  +  4  =  7 

30  (29) 

5 

18 

7 

0 

18 

20 

4  +  3  =  7 

10(10) 

1 

4 

3 

2 

4 

6 

4  +  4  =  8 

4(4) 

2 

1 

1 

0 

2 

2 

Total 

104 

35  (34) 

52  (50) 

14(13) 

3(3) 

71  (68) 

76  (73) 

Note. — Data  are  numbers  of  patients,  with  percentages  in  parentheses.  Gleason  score  ranking  was  performed  with  the  following  algorithm:  3  +  3  =  6  (rank  1),  3  +  4  =  7  (rank  2),  4  +  3  =  7  (rank  3), 
and  4  +  4  =  8  (rank  4). 


staging  at  pathologic  examination  and 
imaging. 

Statistical  Analysis 

The  following  variables  were  record¬ 
ed  for  each  patient:  age,  serum  PSA 
level,  original  surgical  plan,  final  sur¬ 
gical  plan,  biopsy  Gleason  score,  final 
pathologic  Gleason  score,  weeks  be¬ 
tween  biopsy  and  MR  imaging,  weeks 
between  MR  imaging  and  surgery,  and 
components  of  MR  imaging  (spectros¬ 
copy,  diffusion-weighted  imaging,  and 
dynamic  contrast-enhanced  imaging). 
Rinary  variables  were  assigned  sepa¬ 
rately  for  the  right  and  left  side  for 
each  patient  for  both  MR  imaging 
stage  and  surgical  stage  (T2  vs  T3) 
and  for  surgical  margins.  The  type 


of  surgery  performed  (nerve-sparing 
vs  non-nerve-sparing  surgery)  and 
whether  the  planned  surgery  changed 
on  the  basis  of  imaging  findings  was 
also  assigned  a  binary  variable  sepa¬ 
rately  for  the  right  and  left  side.  The 
Gleason  score,  which  was  determined 
at  both  biopsy  and  surgery,  was  an¬ 
alyzed  separately  and  converted  to 
ranked  categories  (see  footnote  in 
Table  2).  Serum  PSA  level  and  patient 
age  were  analyzed  as  continuous  var¬ 
iables.  The  time  between  biopsy  and 
MR  imaging  and  between  MR  imag¬ 
ing  and  surgery  was  also  analyzed  as 
a  continuous  variable.  Rinary  vari¬ 
ables  underwent  multivariate  logistic 
regression  analysis.  Continuous  and 
ranked  variables  underwent  Pearson 


product  moment  and  nonparametric 
Spearman  rank  correlation  analysis 
with  use  of  software  (SPSS  PASW 
17;  SPSS,  Chicago,  Ill).  Although  the 
sample  size  appeared  large  enough 
for  Pearson  product  moment,  because 
we  initially  planned  on  evaluating 
the  Spearman  rank  correlation,  this 
analysis  is  also  presented.  P  <  .05 
was  considered  indicative  of  a  statisti¬ 
cally  significant  difference. 


Results 


Effect  on  Surgery  and  Surgical  Margins 

The  initial  planned  surgical  technique 
was  changed  following  review  of  the 
MR  imaging  data  on  at  least  one  side 
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a.  b. 

Figure  1 :  Example  of  surgical  plan  changed  to  nerve-sparing  technique.  Images  in  59-year-old 
man  with  PSA  level  of  4.3  ng/dL  and  high  volume  at  biopsy  (Gleason  grade  3  +  3  =  6).  Non¬ 
nerve-sparing  surgery  was  initially  planned  on  the  left  on  the  basis  of  biopsy  findings.  MR  imaging 
showed  no  involvement  of  the  neurovascular  bundles  or  seminal  vesicles,  and  surgical  plan  was 
changed  to  bilateral  nerve-sparing  surgery.  (a)T2-weighted  axial  MR  image  shows  low-signal-in- 
tensity  focus  (horizontal  arrow)  at  left  anterolateral  midgland  with  preserved  neurovascular  bundle 
(vertical  arrow),  (b)  Apparent  diffusion  coefficient  map  shows  restricted  diffusion  (arrow),  which 
corresponds  to  the  low-signal-intensity  focus  in  a.  (c)  Dynamic  contrast-enhanced  subtraction 
image  obtained  30  seconds  after  contrast  material  injection  (downsampled  by  a  factor  of  2  in 
both  directions)  shows  early  enhancing  nodule  (arrow),  (d)  Multivoxel  spectroscopic  overlay  shows 
elevated  choline  levels  (arrows)  corresponding  to  anterior  midgland.  Final  pathologic  examination 
demonstrated  Gleason  score  of  3  +  4  =  7  and  confirmed  imaging  findings  of  organ-confined 
cancer,  with  no  suggestion  of  extracapsular  extension  to  the  neurovascular  bundle. 


in  28  of  the  104  patients  (27%)  (Figs 
1-3).  In  17  of  the  28  patients,  sur¬ 
gery  was  changed  from  a  non-nerve¬ 
sparing  technique  to  a  nerve-sparing 
technique — nine  on  the  right  side  only 
and  seven  on  the  left  side  only.  In  one 
patient,  surgery  was  changed  from  a 
non-nerve-sparing  to  a  nerve-sparing 
technique  on  the  right  but  from  a 
nerve-sparing  to  a  non-nerve-sparing 
technique  on  the  left.  In  11  of  the  28 
patients,  surgery  was  changed  from  a 
nerve-sparing  to  a  non-nerve-sparing 
technique — six  on  the  right  side  only, 
four  on  the  left  side  only,  and  one 
bilaterally. 

Surgical  margins  were  positive  in 
five  patients  on  the  right  side  only, 
in  one  patient  on  the  left  side  only, 
and  in  one  patient  bilaterally.  In  two 
of  these  cases,  the  surgical  plan  was 
changed — but  on  the  side  opposite  of 
the  positive  margin.  In  no  case  of  pos¬ 
itive  margins  was  the  surgery  changed 
to  nerve-sparing  surgery  on  the  side 


d. 


of  the  positive  margins.  In  four  addi¬ 
tional  cases,  margins  were  only  focally 
positive  at  the  apex,  which  is  consid¬ 
ered  a  negative  margin  because  this 
has  not  been  shown  to  confer  a  worse 
prognosis  than  organ-confined  disease 
(14)  and  is  managed  by  urologists  as 
a  negative  margin. 

Positive  margins  coincided  with 
pathologic  stage  T3  (pT3)  disease  bilat¬ 
erally  in  the  case  of  bilateral  positive 
margins  and  ipsilaterally  in  one  case 
of  left- sided  positive  margins;  all  other 
cases  of  positive  margins  occurred  in 
stage  pT2  disease,  where  the  surgical 
margin  was  inside  the  prostate  capsule 
and  tumor  was  within  1  mm  of  the 
margin.  Non-nerve -sparing  surgery  was 


C. 


performed  bilaterally  in  the  one  patient 
with  bilateral  positive  margins  and  ip¬ 
silaterally  in  one  patient  with  positive 
margins  on  the  right  side;  in  all  other 
cases  of  positive  margins,  nerve-sparing 
surgery  was  performed  on  the  positive 
side. 

Diffusion-weighted  imaging  was 
performed  in  all  patients  with  posi¬ 
tive  margins  and  dynamic  contrast-en¬ 
hanced  imaging  was  performed  in  five. 
MR  spectroscopy  was  performed  in  all 
patients  except  the  two  with  stage  pT3 
disease. 

Radiologic-Pathologic  Correlation 

In  all  but  one  patient,  tumors  were  either 
posterior  or  both  anterior  and  posterior. 
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c.  d. 

Figure  2:  Example  of  surgical  plan  changed  to  non-nerve-sparing  technique.  Images  in  53-year-old  man 
with  PSA  level  of  1 4  ng/dL  and  biopsy  Gleason  score  of  3  +  3  =  6  in  two  of  six  core  specimens  on  the 
right.  On  the  basis  of  these  findings,  bilateral  nerve-sparing  surgery  was  initially  planned.  After  MR  imaging, 
however,  surgical  plan  was  changed  to  a  non-nerve-sparing  technique  on  the  right,  (a)  Coronal  and  (b)  axial 
T2-weighted  images  show  a  low-signal-intensity  mass  (arrow)  extending  into  the  right  seminal  vesicle,  (c) 
Apparent  diffusion  coefficient  map  shows  restricted  diffusion  (arrow),  corresponding  to  the  low-signal-inten¬ 
sity  focus  on  T2-weighted  images,  (d)  Contrast-enhanced  fat-saturated  T1  -weighted  gradient-echo  image 
shows  enhancing  nodule  (arrow),  corresponding  to  the  area  of  low  signal  intensity  on  T2-weighted  images. 
Final  pathologic  analysis  demonstrated  Gleason  score  of  4  +  3  =  7  and  confirmed  imaging  findings  of  right 
seminal  vesicle  invasion. 


The  one  patient  with  anterior  cancer  had 
stage  T2  disease  at  pathologic  examina¬ 
tion  and  MR  imaging.  Whether  these 
tumors  were  central,  peripheral,  or  both 
was  not  recorded  reliably  by  the  pathol¬ 
ogist;  therefore,  this  was  not  included  in 
the  analysis.  Of  eight  patients  with  stage 
pT3  disease,  disease  was  on  the  left  in 
four  patients,  on  the  right  in  two,  and 
bilateral  in  two.  The  computed  sensitiv¬ 
ity,  specificity,  and  positive  and  negative 
predictive  values  for  predicting  stage  T3 


disease  according  to  side  are  50%  (five 
of  10  cases),  97.5%  (193  of  198  cases), 
50%  (five  of  10  cases),  and  97.5%  (193 
of  198  cases),  respectively  (Table  4).  In 
11  of  the  104  patients,  the  stage  deter¬ 
mined  with  imaging  was  discordant  with 
that  determined  at  pathologic  examina¬ 
tion.  In  one  patient  (who  underwent  both 
diffusion-weighted  and  dynamic  contrast- 
enhanced  imaging),  disease  was  staged 
as  T3b  at  MR  imaging  but  was  T3a  at 
pathologic  examination.  Of  the  remaining 


10  patients,  one  patient  (who  under¬ 
went  MR  spectroscopy  but  not  dynamic 
contrast-enhanced  or  diffusion-weighted 
imaging)  was  classified  with  stage  T3a 
disease  on  the  right  and  stage  T2  disease 
on  the  left  at  MR  imaging  but  was  stage 
pT3a  bilaterally,  four  patients  (one  who 
underwent  MR  spectroscopy  alone,  two 
who  underwent  both  diffusion-weighted 
and  dynamic  contrast-enhanced  imaging, 
and  one  who  underwent  all  three  com¬ 
ponents)  were  classified  with  stage  T2 
disease  at  MR  imaging  and  stage  pT3 
disease  at  pathologic  examination,  and 
the  last  three  patients  (one  who  under¬ 
went  diffusion-weighted  imaging  and  MR 
spectroscopy  and  two  who  underwent 
all  three  components)  were  classified 
with  stage  T3  disease  at  MR  imaging  but 
stage  pT2  disease  at  pathologic  exami¬ 
nation.  Of  the  mismatches,  two  patients 
had  positive  margins — one  patient  was 
understaged  as  having  stage  T2  disease 
with  microscopic  extracap sular  extension 
at  MR  imaging  (diffusion-weighted  and 
dynamic  contrast-enhanced  imaging)  and 
one  was  determined  to  have  stage  T2 
disease  on  the  side  of  the  positive  mar¬ 
gin  at  both  MR  imaging  and  pathologic 
examination  (discordant  contralaterally) . 
It  should  be  noted  that  “possible  micro¬ 
scopic  extracapsular  extension”  at  MR 
imaging  was  recorded  for  47  patients. 
Finally,  of  173  sides  in  which  the  origi¬ 
nal  plan  was  for  nerve-sparing  surgery, 
168  (97.1%)  were  classified  as  stage  T2 
at  MR  imaging.  Of  these  168  sides,  163 
(97%)  were  pathologically  proved  to  be 
stage  pT2.  The  breakdown  of  patients  ac¬ 
cording  to  surgical  plan  and  surgical  mar¬ 
gin  status  is  provided  in  Figure  4. 

Statistical  Analysis 

Analysis  of  imaging  stage  according 
to  side  found  a  significant  association 
with  surgical  stage  and  a  change  in 
surgery  bilaterally  but  a  significant 
association  with  stage  pT3  disease 
and  positive  margins  on  the  left  only 
(Table  5).  Logistic  regression  analysis 
of  positive  margins  showed  no  asso¬ 
ciation  with  change  in  surgery  and 
a  significant  relationship  with  imag¬ 
ing  and  surgical  stage  on  the  left  but 
not  on  the  right.  No  association  was 
found  between  absence  of  any  of  the 
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Figure  3:  Example  of  surgical  plan  changed  to  nerve-sparing  technique.  Images  in 
67-year-old  man  with  PSA  level  of  5.5  ng/dL  and  biopsy  Gleason  score  of  3  +  4  =  7  in 
seven  of  eight  core  specimens.  Patient  was  originally  scheduled  for  non-nerve-sparing 
surgery  on  the  left,  but  after  MR  imaging  the  surgical  plan  was  changed  to  nerve-sparing 
surgery,  (a)  Axial  T2-weighted  image  shows  diffuse  low  signal  intensity  bilaterally  but 
an  intact  capsule  (arrows),  (b)  Axial  apparent  diffusion  coefficient  map  shows  restricted 
diffusion  (arrow),  corresponding  to  the  left  peripheral  gland,  that  was  also  confined  to  the 
organ,  (c)  Axial  contrast-enhanced  fat-saturated  T1  -weighted  gradient-echo  image  shows 
enhancing  nodule  (arrow)  corresponding  to  areas  of  restricted  diffusion,  again  confined 
to  the  organ,  (d)  Multivoxel  spectroscopic  overlay  shows  elevated  choline  level  (arrows) 
corresponding  to  the  left  (also  right)  rectoprostatic  angle  midgland.  Final  pathologic 
examination  demonstrated  Gleason  score  of  3  +  4  =  7  and  confirmed  imaging  findings  of 
organ-confined  cancer  with  no  suggestion  of  extracapsular  extension  to  the  neurovascu¬ 
lar  bundle. 


components  of  MR  imaging  (diffusion- 
weighted  imaging,  dynamic  contrast- 
enhanced  imaging,  or  MR  spectros¬ 
copy)  and  margins,  grade,  stage  (at 
imaging  or  pathologic  examination), 
surgery  type,  or  change  in  surgical 
plan  owing  to  imaging  findings. 

A  significant  Pearson  correlation 
was  found  between  surgical  grade  and 
both  serum  PSA  level  and  stage  pT3 
disease  and  between  biopsy  grade  and 
serum  PSA  level  but  not  stage  pT3 
disease,  and  neither  had  a  significant 
correlation  with  age  or  change  in  sur¬ 
gery  (Table  6).  A  significant  correla¬ 
tion  was  also  found  between  age  and 
change  in  surgery.  At  Spearman  rank 
coefficient  analysis,  we  found  a  signifi¬ 
cant  correlation  between  surgical  grade 
and  PSA  level  and  stage  pT3  disease 
as  well  as  change  in  surgery  and  be¬ 
tween  biopsy  grade  and  PSA  level  but 
none  of  the  other  parameters.  Again, 


a  significant  correlation  was  found  be¬ 
tween  age  and  change  in  surgery. 


Discussion 


Approximately  40%  of  patients  with 
localized  prostate  cancer  choose  some 
form  of  surgical  resection  for  treatment 

(15) .  In  2008,  approximately  60000 
patients  opted  to  undergo  RALP  rather 
than  conventional  open  prostatectomy 

(16) .  In  any  surgical  approach,  surgeons 
must  balance  the  desire  to  achieve  a 
cancer-free  resection  margin  with  the 
need  to  minimize  postoperative  mor¬ 
bidity,  which  with  prostate  cancer  sur¬ 
gery  includes  incontinence,  impotence, 
and  erectile  dysfunction.  The  neurovas¬ 
cular  bundles,  which  mediate  erectile 
function  and  continence  (17),  lay  pos¬ 
terolateral  to  the  prostatic  capsule  and 
adjacent  to  the  peripheral  zone,  where 


70%  of  cancers  arise.  In  open  radical 
prostatectomy,  surgeons  typically  iden¬ 
tify  and  spare  both  neurovascular  bun¬ 
dles  to  preserve  patient  potency  (18). 
However,  a  surgeon  who  believes  that 
an  area  beyond  the  capsule  is  suspi¬ 
cious  for  cancer  can  widely  resect  the 
neurovascular  bundle  and  surrounding 
tissue  to  achieve  negative  surgical  mar¬ 
gins.  By  resecting  the  neurovascular 
bundle,  the  risk  of  postoperative  impo¬ 
tence  increases  substantially  (19). 

During  RALP,  surgeons  lack  the  tac¬ 
tile  feedback  they  rely  on  during  open 
surgery  to  determine  the  extent  of  ma¬ 
lignant  involvement  and  resection,  po¬ 
tentially  affecting  surgical  outcome.  As 
in  open  radical  prostatectomy,  surgeons 
performing  RALP  routinely  dissect  the 
neurovascular  bundles  off  the  pro  static 
capsule  and  spare  the  nerves — but 
without  haptic  feedback.  Unlike  clinical 
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Figure  4 


Figure  4:  Diagram  shows  the  relative  number  of  patients  stratified  according  to  surgery,  imaging  findings, 
and  surgical  margins.  neg=  Negative,  pos=  positive. 


Table  4 


Imaging  versus  Pathologic  Stage 

Pathologic  Stage 


Parameter*  T2  T3 


MR  imaging  stage 

T2 

193 

5 

T3 

5 

5 

Sensitivity  (%) 

50 

Specificity  (%) 

97.5 

PPV  (%) 

50 

NPV  (%) 

97.5 

*  NPV  =  negative  predictive  value,  PPV  =  positive 
predictive  value. 


variables  (PSA  level,  findings  from  dig¬ 
ital  rectal  examination),  results  from 
prostate  MR  imaging  are  spatially  local¬ 
ized  and  allow  surgeons  to  individually 
sculpt  the  extent  of  surgical  resection 
(20).  Prostate  MR  imaging  data  map 
an  individual  cancer’s  relationship  to 
the  capsule  and  neurovascular  bundles, 
thus  assisting  in  the  determination  of 
performing  nerve-sparing  versus  non¬ 
nerve-sparing  surgery.  A  priori  knowl¬ 
edge  of  extracapsular  extension  would 
allow  the  surgeon  to  plan  non-nerve  - 
sparing  resection  of  the  neurovascu¬ 
lar  bundle  on  the  involved  side  (21). 
Conversely,  sparing  the  neurovascular 
bundles  when  extracapsular  extension 
is  present  increases  the  probability  of 
a  positive  surgical  margin,  the  need  for 
postoperative  radiation  therapy,  and  lo¬ 
cal  cancer  recurrence  (22). 

In  patients  undergoing  open  rad¬ 
ical  prostatectomy,  MR  imaging  has 
been  shown  to  improve  the  accuracy 
of  the  surgeon’s  decision  to  resect  or 
preserve  the  neurovascular  bundles 
(20) .  Our  study  differs  from  this  land¬ 
mark  study  in  several  ways.  First,  our 
study  could  be  viewed  by  some  as  be¬ 
ing  more  consistent  with  true  clinical 
practice  because  the  surgical  plan  was 
informed  by  the  report  and  online  im¬ 
ages  but  rarely  by  direct  radiologist 
consultation.  In  addition,  not  all  pa¬ 
tients  underwent  MR  spectroscopy. 
The  main  reason  why  MR  spectros¬ 
copy  was  not  performed  in  all  patients 
is  because  it  is  not  covered  by  most 
insurance  carriers  and  not  all  patients 


can  pay  for  the  additional  cost.  An¬ 
other  difference  between  studies  is 
that  our  scoring  of  stage  T3  disease  at 
MR  imaging  and  the  surgical  plan  was 
binary,  not  quinary.  Although  this  pre¬ 
cludes  receiver  operating  characteris¬ 
tic  analysis,  the  surgeon’s  decision  to 
spare  or  resect  neurovascular  bundles 
was,  in  the  context  of  this  analysis, 
binary;  therefore,  it  was  requested 
that  our  assessment  of  stage  at  imag¬ 
ing  reflect  this.  In  addition,  all  of  the 
images  in  our  study  were  read  by  the 
same  two  radiologists  to  reduce  the 
likelihood  of  interobserver  bias;  in  the 
previous  study,  all  images  were  read 
by  one  of  two  radiologists.  Last,  an 
important  point  is  that  we  looked  at 
how  MR  imaging  changed  the  surgical 
plan  and  the  final  pathologic  diagno¬ 
sis,  not  just  the  presence  of  positive 
surgical  margins  as  in  the  previous 
study.  The  correlation  of  imaging  and 
pathologic  findings  was  not  our  pri¬ 
mary  end  point;  rather,  we  wanted 
to  evaluate  the  impact  of  MR  imaging 


data  on  the  extent  of  surgical  resec¬ 
tion  bilaterally.  This  last  aspect  is  the 
most  important  because  the  tactile 
response  of  the  surgeon  using  robotic 
technology  is  limited,  increasing  the 
reliance  on  imaging  for  determin¬ 
ing  stage.  This  informs  the  selection 
of  patients  because  this  technique  is 
optimized  for  preservation  of  neuro¬ 
vascular  bundles  and  is  therefore  of 
less  value  in  gross  T3  disease  (16).  In 
the  previous  study  (20),  all  patients 
underwent  open  prostatectomy  and 
all  underwent  identical  imaging  with¬ 
out  diffusion-weighted  or  dynamic 
contrast-enhanced  pulse  sequences. 
Analysis  in  this  previous  study  used 
histopathologic  correlation  of  neuro¬ 
vascular  bundle  invasion  as  a  primary 
end  point  for  modeling  risk  stratifica¬ 
tion  rather  than  analyzing  the  effect 
on  surgical  technique  and  margins.  In 
addition,  the  surgical  plan  was  altered 
in  39%  of  the  neurovascular  bundles 
in  that  study,  which  is  higher  than 
the  percentage  we  found.  This  may 
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Table  6 


Pearson  Product  Moment  and  Spearman  Rank  Correlation  Coefficients 

Parameter 

Change  in  Surgery 

Surgical  Grade 

Biopsy  Grade 

Stage  T3 

PSA  Level 

Pearson  analysis 

Age 

0.263  (.007) 

0.078  (.429) 

0.007  (.941) 

0.074  (.453) 

0.037  (.712) 

PSA  level 

0.007  (.941) 

0.295  (.003) 

0.217  (.028) 

0.141  (.156) 

Stage  T3  disease 

0.013  (.900) 

0.247  (.011) 

0.143  (.147) 

Biopsy  grade 

0.141  (.153) 

Surgical  grade 

0.174  (.077) 

Spearman  analysis 

Age 

0.234  (.017) 

0.008  (.937) 

0.012  (.904) 

0.073  (.459) 

0.056  (.577) 

PSA  level 

0.029  (.772) 

0.334  (.001) 

0.213  (.031) 

0.234  (.018) 

Stage  T3  disease 

0.013  (.900) 

0.222  (.023) 

0.130  (.187) 

Biopsy  grade 

0.121  (.222) 

Surgical  grade 

0.199  (.043) 

Note. — Numbers  in  parentheses  are  P  values. 


be  due  in  part  to  patient  selection 
and  the  emphasis  on  nerve-sparing 
technique  placed  on  robotic-assisted 
surgery.  Of  173  sides  in  which  the 
original  plan  was  for  nerve-sparing 
surgery,  168  (97.1%)  were  stage  T2 
at  MR  imaging.  Of  these  168  sides, 
163  (97%)  were  proved  to  be  stage 
pT2  at  pathologic  examination. 

In  28  of  our  104  patients  (27%),  the 
surgical  plan  was  changed  because  of  in¬ 
formation  obtained  only  at  prostate  MR 
imaging;  17  of  those  28  patients  (61%) 
underwent  nerve-sparing  surgery.  On 
the  basis  of  the  results  of  preoperative 
MR  imaging  and  MR  spectroscopy,  pa¬ 
tients  were  able  to  undergo  a  less  mor¬ 
bid  nerve-sparing  surgery  without  com¬ 
promising  oncologic  outcome.  Although 
the  overall  rate  of  positive  surgical  mar¬ 
gins  was  6.7%  (seven  of  104  patients), 
86%  (six  of  seven  patients)  of  these 
patients  had  undergone  a  more  aggres¬ 
sive  non-nerve-sparing  technique  with 
concordant  MR  findings  on  the  affected 
side.  In  four  additional  patients,  posi¬ 
tive  margins  were  only  focally  positive 
at  the  apex;  this  has  not  been  shown  to 
confer  a  worse  prognosis  than  organ- 
confined  disease  (14). 

The  specificity  of  prostate  MR 
imaging  in  the  differentiation  of  T2 
disease  from  T3  disease  was  97.5%. 
This  suggests  that,  in  our  patient  pop¬ 
ulation,  if  a  patient’s  prostate  MR  im¬ 
age  shows  T2  disease,  a  nerve-sparing 


technique  can  be  used  with  a  high 
degree  of  safety  and  without  compris¬ 
ing  oncologic  outcome.  Our  sensitivity 
and  specificity  were  higher  than  that 
reported  in  a  recent  retrospective 
analysis  of  endorectal  coil  MR  imag¬ 
ing  by  Rrajtbord  et  al  (23)  (50%  and 
97.5%  vs  43%  and  73%,  respectively). 
However,  Rloch  et  al  (24)  showed  a 
similar  specificity  of  95%. 

Our  study  had  several  limitations. 
First,  our  final  outcome  “standard  of 
reference”  was  positive  surgical  mar¬ 
gins,  not  clinical  outcome.  Follow-up 
studies  looking  specifically  at  these  var¬ 
iables  are  required  to  determine  how 
the  refinement  of  surgical  technique 
guided  by  MR  imaging  findings  affects 
patient  continence  and/or  potency  and 
cancer  recurrence  rates.  Second,  the 
final  pathology  reference  was  step- 
section  histologic  examination  and  not 
whole-mount  histologic  examination 
because  this  was  not  performed  at  our 
institution  during  this  time.  Third,  our 
results  represent  one  surgeon’s  expe¬ 
rience  with  a  consecutive,  nonselected 
patient  population  and  may  not  be  rep¬ 
resentative  of  the  experience  at  other 
centers.  Fourth,  our  studies  were  read 
in  consensus  and  therefore  do  not  re¬ 
flect  clinical  practice.  Also,  most  clin¬ 
ical  and  research  studies  are  read  by 
one  of  a  core  group  of  radiologists. 
Fifth,  because  of  our  study  design 
(provide  the  highest  level  of  imaging 


for  optimal  patient  care,  such  that  all 
patients  underwent  diffusion-weighted 
and  dynamic  contrast-enhanced  imag¬ 
ing  but  only  once  available),  we  could 
not  control  for  protocol.  Finally,  the 
protocol  of  the  MR  examinations  var¬ 
ied,  which  precluded  evaluation  of  in¬ 
dividual  components  of  MR  imaging. 
Despite  these  limitations,  our  compar¬ 
ison  of  MR  imaging  to  the  standard  of 
reference  of  histopathologic  examina¬ 
tion  yielded  excellent  results.  Possibly 
the  most  surprising  finding  was  that 
not  only  were  no  cases  of  positive  mar¬ 
gins  found  on  the  side  where  MR  imag¬ 
ing  resulted  in  a  change  to  the  surgical 
plan,  but  in  only  one  of  seven  cases 
where  imaging  found  stage  T3  disease 
was  the  surgical  plan  changed.  This 
may  suggest  that  imaging  did  not  add 
to  the  decision  because  non-nerve¬ 
sparing  surgery  was  already  planned 
in  six  of  these  cases.  Overall,  however, 
imaging  resulted  in  a  change  to  the 
surgical  plan  in  28  of  the  104  patients 
(27%),  with  most  surgeries  changed  to 
a  nerve-sparing  technique;  in  addition, 
there  were  no  positive  margins  on 
the  side  where  the  surgical  plan  was 
changed  to  a  nerve-sparing  technique, 
potentially  decreasing  the  chance  of 
postsurgical  morbidity. 

To  our  knowledge,  this  is  the  first 
report  demonstrating  the  value  of  MR 
imaging  in  planning  RALP.  On  the  basis 
of  our  initial  experience,  prostate  MR 
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imaging  may  be  useful  for  helping  sur¬ 
geons  plan  the  extent  and  side  of  nerve 
sparing  during  RALP. 
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Application  of  Compressed  Sensing  to  Multidimensional 
Spectroscopic  Imaging  in  Human  Prostate 

Jon  K.  Furuyama,  Neil  E.  Wilson,  Brian  L.  Burns,  Rajakumar  Nagarajan,  Daniel  J.  Margolis, 
and  M.  Albert  Thomas* 


The  application  of  compressed  sensing  is  demonstrated  in 
a  recently  implemented  four-dimensional  echo-planar  based 
J-resolved  spectroscopic  imaging  sequence  combining  two 
spatial  and  two  spectral  dimensions.  The  echo-planar  readout 
simultaneously  acquires  one  spectral  and  one  spatial  dimen¬ 
sion.  Therefore,  the  compressed  sensing  undersampling  is 
performed  along  the  indirectly  acquired  spatial  and  spectral 
dimensions,  and  the  reconstruction  is  performed  using  the 
split  Bregman  algorithm,  an  efficient  TV-minimization  solver. 
The  four-dimensional  echo-planar-based  J-resolved  spectro¬ 
scopic  imaging  data  acquired  in  a  prostate  phantom  contain¬ 
ing  metabolites  at  physiological  concentrations  are  accurately 
reconstructed  with  as  little  as  20%  of  the  original  data.  Experi¬ 
mental  data  acquired  in  six  healthy  prostates  using  the  external 
body  matrix  “receive”  coil  on  a  3T  magnetic  resonance  imag¬ 
ing  scanner  are  reconstructed  with  acquisitions  using  only  25% 
of  the  Nyquist-Shannon  required  amount  of  data,  indicating  the 
potential  for  a  4-fold  acceleration  factor  in  vivo,  bringing  the 
required  scan  time  for  multidimensional  magnetic  resonance 
spectroscopic  imaging  within  clinical  feasibility.  Magn  Reson 
Med  67:1499-1505,  2012.  ©  2012  Wiley  Periodicals,  Inc. 

Key  words:  Compressed  sensing;  prostate;  magnetic  reso¬ 
nance;  citrate;  echo-planar  J-resolved  spectroscopic  imaging 

Magnetic  resonance  spectroscopy  and  magnetic  resonance 
spectroscopic  imaging  have  evolved  as  powerful  research 
tools  for  their  ability  to  study  the  underlying  biochem¬ 
istry  of  tissue  (1)  and  can  greatly  complement  standard 
magnetic  resonance  imaging  (MRI).  One  of  the  major  draw¬ 
backs  of  one-dimensional  (ID)  spectroscopy  is  the  inherent 
overcrowding  of  spectra  due  to  overlapping  peaks.  This 
limitation  can  be  alleviated  with  the  addition  of  more 
spectral  dimensions  by  which  resonances  can  be  spread 
apart,  increasing  the  spectral  dispersion  (2).  Different  two- 
dimensional  (2D)  spectroscopic  techniques  have  been  suc¬ 
cessfully  used  in  vivo  such  as  /-resolved  spectroscopy 
(3)  and  localized  correlated  spectroscopy  (4).  However, 
these  were  originally  limited  to  single  voxel  acquisitions. 
To  increase  the  spatial  coverage,  the  localized  correlated 
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spectroscopy  and  J-resolved  spectroscopy  sequences  were 
recently  modified  with  an  echo-planar  spectroscopic  imag¬ 
ing  (EPSI)  (5-7)  readout  to  yield  2D  spectra  from  multiple 
voxels  in  a  single  experiment,  called  echo-planar  corre¬ 
lated  spectroscopic  imaging  (8)  and  echo-planar  /-resolved 
spectroscopic  imaging  (EP-JRESI)  (9),  respectively.  Despite 
the  rapid  acquisition  of  EPSI,  such  four-dimensional  (4D) 
scans  still  require  a  considerable  amount  of  scan  time  (~20 
40  mins/average  depending  on  the  desired  spatial/spectral 
resolution),  severely  limiting  clinical  applicability. 

In  recent  years,  the  field  of  compressed  sensing  (CS) 
has  garnered  much  interest  in  the  imaging  community  for 
its  ability  to  reconstruct  images  from  datasets  whose  sam¬ 
pling  does  not  meet  the  Nyquist-Shannon  criterion  (10,11). 
CS  operates  under  the  assumption  that  the  fully  sampled 
data  is  sparse  within  some  transform  domain.  The  CS 
reconstruction  attempts  to  enforce  this  assumption  in  that 
particular  transform  domain,  while  maintaining  fidelity 
with  the  acquired  measurements.  Since  the  application 
of  CS  in  MRI  was  demonstrated  (12,13),  there  has  been  a 
trove  of  proposed  applications  for  imaging  studies  as  well 
as  high  resolution  nuclear  magnetic  resonance,  many  of 
which  hold  the  potential  to  reduce  scan  times  by  a  factor 
of  2-8  times  (14-18). 

Despite  the  many  potential  applications  in  standard  MRI, 
there  have  not  been  many  proposed  implementations  for 
use  in  magnetic  resonance  spectroscopic  imaging  (19).  Part 
of  the  reason  is  that  while  CS  can  provide  acceleration 
factors  ranging  from  2  to  8  times  in  standard  chemical 
shift  imaging  with  multiple  phase  encoding  directions,  the 
acceleration  factor  of  using  an  EPSI  readout  is  on  the  order 
of  16-32.  The  challenge  is  thus  to  further  accelerate  EPSI 
with  CS.  Since  EPSI  simultaneously  reads  out  one  spatial 
and  one  spectral  dimension,  there  remains  only  one  dimen¬ 
sion  that  can  be  undersampled  in  a  2D  spatial  acquisition. 
However,  it  has  been  previously  discussed  that  CS  can  per¬ 
form  better  when  the  undersampling  is  spread  across  mul¬ 
tiple  dimensions  (20).  Hu  et  al.  employed  pseudo-random 
phase-encoding  blips  during  the  EPSI  readout  to  create 
nonuniform  sampling  along  the  spatial  as  well  as  spec¬ 
tral  dimensions  for  use  in  hyperpolarized  13C  spectroscopic 
imaging  (19). 

In  the  EP-JRESI  sequence,  the  EPSI  readout  simultane¬ 
ously  acquires  one  spatially  encoded  dimension  ( kx )  and 
one  temporal  dimension  [t2),  leaving  the  remaining  spa¬ 
tial  and  spectral  dimensions  ( ky  and  fi,  respectively)  to 
be  incrementally  collected.  We  propose  the  use  of  nonuni¬ 
form  undersampling  (NUS)  in  the  remaining  kyt\  plane, 
using  CS  to  reconstruct  the  equivalent  missing  data  to  a 
fully  sampled  4D  EP-JRESI  acquisition.  Despite  the  mixing 
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of  spatial  and  spectral  dimensions  in  the  reconstruction, 
the  sparsity  requirement  for  reconstruction  is  shown  to 
still  be  satisfied,  as  required  by  CS.  Since  the  kyt 1  plane 
is  incrementally  acquired,  NUS  can  be  trivially  applied  to 
the  sequence  without  pseudo-random  gradient  blips. 

Using  fully  sampled  prostate  phantom  datasets,  the  NUS 
can  be  simulated  and  the  data  can  be  reconstructed  using 
publicly  available  U-norm  minimization  algorithms.  We 
show  that  the  CS  reconstruction  performs  favorably  with  as 
little  as  20-25%  of  the  original  data,  implying  an  accelera¬ 
tion  factor  of  4-5  times  and  scan  times  less  than  4-5  mins 
for  the  NUS  4D  EP-JRESI  sequence.  Nonuniformly  under¬ 
sampled  data  were  collected  in  the  prostate  of  healthy 
volunteers  amounting  to  25%  of  the  sampling  required 
by  the  Nyquist-Shannon  criterion.  The  reconstructed  NUS 
data  shows  good  spatial  and  spectral  quality  in  the  under¬ 
sampled  dimensions,  comparable  to  the  fully  sampled 
reconstruction. 

THEORY 

In  order  for  CS  to  be  successfully  applied,  the  follow¬ 
ing  criteria  must  be  satisfied.  First,  the  data  must  have 
a  sparse  representation:  namely,  there  must  exist  some 
transform  domain  in  which  the  data  has  many  coefficients 
that  are  zero  or  nearly  zero.  Second,  the  NUS  must  be 
such  that  it  produces  incoherent  aliasing  artifacts  within 
the  transform  domain,  which  is  easily  obtained  by  sam¬ 
pling  the  dataset  at  random.  With  these  conditions,  the 
fully  sampled  dataset  can  be  recovered  using  nonlinear 
reconstruction  algorithms  that  attempt  to  enforce  the  spar¬ 
sity  while  simultaneously  maintaining  the  fidelity  of  the 
original  measurements  to  within  the  noise.  CS  solves  the 
constrained  optimization  problem  (12) 

min  ||'I'm||1  s.t.||.Fuin  -  y\\%  <  €2  [1] 

m 

where  is  the  transform  operator  in  which  the  recon¬ 
structed  data,  m,  is  sparse,  Tu  is  the  undersampled  Fourier 
transform,  y  is  the  sampled  data,  e  is  a  fidelity  factor,  and 
||x|U  is  the  in  norm:  Hxll^  =  (JA  \xi\n)1/n,  for  n  >  0.  The 
minimization  of  the  fh-norm  in  Eq.  1  promotes  sparsity 
by  trying  to  reduce  the  number  of  nonzero  coefficients  as 
much  as  possible.  The  U-norm  objective  function  is  offset 
by  the  l2-noiin  constraint  that  the  reconstruction  remains 
consistent  with  the  measured  data  within  the  tolerance,  e. 

In  MRI,  wavelets  have  been  the  popular  transform 
domain  in  which  the  U-norm  is  minimized.  However,  in 
this  work,  the  data  is  undersampled  in  the  kyt\  plane  in 
which  the  wavelet  transform  may  not  necessarily  be  as 
optimal  as  in  MRI.  As  described  by  Lustig  et  al.,  mini¬ 
mizing  finite  differences  (4>  =  V),  or  total  variation  (TV), 
can  also  be  used  alongside  the  wavelet  transform  as  a 
form  of  smoothing  the  data  (12).  Other  work  has  shown 
that  minimizing  TV  alone  is  sufficient  to  adequately  recon¬ 
struct  undersampled  datasets  (13,21,22).  Here,  we  use  TV 
exclusively  to  sparsify  the  data. 

The  constrained  problem  in  Eq.  1  can  alternatively  be 
written  as  an  unconstrained  problem 

X  ? 

nunllVy^mlla  +  -\\Fum-y  ||* 


where  X  is  a  regularization  parameter  that  weighs  the  spar¬ 
sity  against  the  data  consistency,  m  =  S[x,y,F1,F2 )  is  the 
final  dataset,  Tu  is  applied  only  along  the  y  and  F1  dimen¬ 
sions,  y  =  s[x,ky,t1,F2)  is  the  sampled  data,  and  Vyya  is 
the  gradient  in  the  yF1  plane. 

In  a  fully  sampled  EP-JRESI  sequence,  the  final  4D 
dataset,  S[x,y,F1,F2)  =  jr{s(kx,k7,  t1}  t2)},  is  composed  of 
two  spatial  and  two  spectral  dimensions.  The  simplest  way 
to  visualize  the  data  is  to  consider  that  for  every  point  in 
the  xy  spatial  plane  (voxel),  there  is  a  2D  /-resolved  spec¬ 
trum  [F1,F2).  An  example  of  this  is  shown  in  Fig.  la,b. 
Experimentally,  the  data  is  collected  such  that  the  entire 
kxt2  plane  is  acquired  during  a  single  EPSI  readout  for  a 
given  ( ky ,  fa)  point.  Nonuniformly  undersampling  the  kyt 1 
plane  will  result  in  incoherent  aliasing  in  the  yF1  plane  fol¬ 
lowing  Fourier  transform.  Figure  lc  shows  the  yF1  plane 
for  the  choline  peak  at  3.2  ppm  (dotted  line  in  IB)  for 
a  given  point  along  the  x-direction  (dotted  line  in  1A). 
Figure  lc  also  shows  the  finite  differences  along  the  y  and 
F1  dimensions.  While  it  can  be  seen  that  the  yF1  plane 
appears  sparse  within  the  identity  domain  (4>  =  1),  it  can 
be  seen  that  || Vm||i  <  ||lm||i,  as  sought  out  by  Eq.  1,  and 
thus  minimizing  the  finite  differences  are  well  suited  for 
CS  reconstruction.  The  optimal  sparsifying  transform  for 
a  given  application  is  not  yet  known  and  is  the  subject  of 
ongoing  research  (23). 

METHODS 

Phantom 

To  determine  the  feasibility  as  well  as  evaluate  the  per¬ 
formance  of  the  CS  reconstruction,  numerous  retrospec¬ 
tive  undersampling  and  reconstructions  were  performed 
on  a  fully  sampled  prostate  EP-JRESI  phantom  dataset. 
A  500  mL  prostate  phantom  was  prepared  containing  the 
following  metabolites  at  physiological  concentrations  as 
reported  in  healthy  human  prostate  (24,25):  citrate  (Cit, 
50  mM),  creatine  (Cr,  5mM),  choline  (Cho,  ImM),  sper¬ 
mine  (Spm,  6mM),  myo-inositol  (ml,  10  mM),  phospho- 
choline  (PCh,  2mM),  taurine  (Tau,  3mM),  glutamate 
(Glu,  4mM),  glutamine  (Gin,  2.5  mM)  and  scyllo-inositol 
(Scy,  0.8  mM).  The  dataset  was  localized  with  a  field  of 
view  (FOV)  of  16  x  16  cm2  on  a  16  x  16  grid  with  a 
slice  thickness  of  2  cm  for  an  individual  voxel  volume 
of  2  cm3;  512  bipolar  gradient  echo  pairs  ( t2  samples) 
were  collected  with  a  repeat  time  of  0.84  ms  and  100fi 
increments  with  Afa  =  1ms  resulting  in  F1  and  F2 
bandwidths  of  1000  Hz  and  1190  Hz,  respectively.  With 
pulse  repetition  time/echo  time  =  1500/30  ms  and  one 
average,  the  EP-JRESI  phantom  scan  duration  was  40  mins. 
Water  suppression  was  performed  using  the  WET  sequence 
(26)  just  before  PRESS  localization  (27).  To  correct  for 
eddy-currents  generated  by  the  EPSI  readout,  a  nonwater- 
suppressed  scan  was  acquired  as  a  reference  (28). 

The  NUS  was  simulated  by  zeroing  data  points  in  the 
fully  sampled  kyt\  plane  according  to  the  exponentially 
decaying  sampling  density 

,  a  f  \ky\  t\ 

p{ky,t i)  =  exp  -  — 


[2] 


[3] 
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FIG.  1.  a:  71 -weighted  axial  MRI  of  a  prostate  phantom  where  the 
white  box  indicates  the  PRESS  localized  volume,  b:  2D  J-resolved 
spectrum  from  the  center  voxel  (yellow  box  in  a)  of  the  prostate  phan¬ 
tom,  c:  Demonstration  of  the  sparsity  of  the  undersampled  yF-i  plane, 
S(y,F-|),  for  a  given  x  and  F2  point  as  indicated  by  the  dotted  lines 
in  (a)  and  (b)  as  well  as  the  finite  differences  of  the  same  plane, 
vy,FS(y,F  i). 


where  p [ky,ti)  is  the  probability  a  data  point  is  sampled, 
—ky> max  <  ky  <  ky>max,  0  <  t\  <  and  a  and  b  are  flex¬ 

ible  parameters  that  determine  the  percentage  of  data  that 
is  sampled.  The  signal  in  the  indirect  t\  dimension  has  an 
exponential  decay  envelope  determined  by  T2 -weighting 
and  the  the  center  of  k- space  has  the  highest  signal  inten¬ 
sity.  Therefore,  Eq.  3  ensures  that  the  regions  in  the  kyt\ 
plane  with  the  greatest  signal-to-noise  ratio  (SNR)  are  col¬ 
lected.  Different  trials  were  conducted  for  NUS,  keeping 
50%,  33%,  25%,  20%,  16%,  and  10%  of  the  original  phan¬ 
tom  data,  corresponding  to  acceleration  factors  from  2  to  10 
times.  For  each  simulation,  a  different  mask  was  randomly 
generated  with  a  =  b  such  that  the  degree  of  undersampling 
was  distributed  equally  between  the  ky  and  t\  dimensions. 

In  Vivo 

The  performance  of  the  reconstruction  algorithm  was 
tested  in  vivo  by  collecting  undersampled  data  in  the 
prostate  of  6  healthy  volunteers  (age  range  25-57  years  old) 
using  a  body  matrix  “receive  assembly”  on  a  Siemens  3T 
Trio-TIM  MRI  scanner  (Siemens  Medical  Solutions,  Erlan¬ 
gen,  Germany)  running  the  VBl7a  platform.  The  FOV  was 
a  16  x  16  cm2  area  localized  onto  a  16  x  16  grid  with  slice 
thickness  of  2  cm,  for  an  individual  voxel  volume  of  2  cm3. 
The  usual  fully  sampled  data  consists  of  64  t\  increments. 
The  scanner  was  programmed  to  collect  only  25%  of  the 
fully  sampled  data  matrix  according  to  the  sampling  den¬ 
sity  in  Eq.  3.  Points  that  were  not  collected  were  set  to  zero. 
Two  averages  were  acquired  for  a  total  scan  duration  of 
12  min  48  s.  A  nonwater-suppressed  scan  was  also  acquired 
as  a  reference  scan  for  eddy-current  correction  and  coil 
combination  purposes  by  fully  sampling  ky  with  only  the 
first  fi  increment,  adding  30  s  to  the  total  scan  duration. 

Reconstruction 

All  data  reconstructions  were  performed  by  solving  Eq.  2 
using  the  split  Bregman  iterative  method  (29)  with  X  =  1 
such  that  the  data  fidelity  is  weighed  equally  against  the 
minimization  of  the  finite  differences.  The  algorithm  was 
iterated  until  V7;^1  m  was  minimized  while  maintaining 
a  normalized  error  of  no  greater  than  1  part  per  million 
(ppm),  namely  ||  JFum  —  y  |||  <  lx  10-6.  While  the  recon¬ 
struction  was  performed  on  the  entire  dataset,  the  TV 
operator,  Vy^,  only  acted  on  data  in  the  undersampled 
kyt\  plane.  The  entire  4D  reconstruction  required  roughly 
30mins/channel  on  a  standard  desktop  PC. 

For  both  phantom  and  in  vivo  scans,  we  have  used  the 
maximum-echo  sampling  method  as  reported  by  Schulte 
et  al.  (30).  This  yields  a  spectrum  with  a  COSY-like  diagonal 
where  the  /-resolved  peaks  are  centered  around  the  diag¬ 
onal.  Each  t\  row  was  time-shifted  during  postprocessing 
such  that  the  diagonal  is  rotated  45°,  having  the  appearance 
of  a  conventional  2D  /-resolved  spectrum  (31)  with  an  F1 
bandwidth  of  ±250  Hz. 

For  the  phantom,  the  CS  reconstruction  was  applied  to 
a  processed  fully-sampled  dataset  that  had  NUS  imposed 
retrospectively.  However,  for  the  in  vivo  scan  more  steps 
were  required  since  the  raw  metabolite  data  itself  was 
undersampled.  A  nonwater-suppressed  scan  was  collected 
to  correct  for  eddy  current  distortions  and  to  serve  as 
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FIG.  2.  a:  Spatial  distribution  of  the  fully 
sampled  J-resolved  citrate  multiplet  at 
2.6  ppm  on  the  7"i  -weighted  axial  MRI. 
b:  The  effect  of  the  incoherent  aliasing 
along  both  the  y-direction  as  well  as 
the  F-i -dimension  as  a  function  of  the 
percentage  of  the  nonuniformly  sampled 
data.  For  each  amount  of  kept  data,  the 
column  along  the  y-direction  for  a  given 
point  along  x  (left  column)  is  marked  by 
the  dotted  yellow  box  in  (a),  and  the  F\ 
column  (right  column)  is  centered  around 
the  citrate  multiplet  at  F2  =  2.6  ppm. 


a  reference  for  coil  combination.  The  4D  reconstruction 
was  performed  on  the  raw  metabolite  data  for  each  chan¬ 
nel  individually,  followed  by  eddy-current  correction  and 
coil  combination.  For  a  quantitative  comparison  of  the  in 
vivo  data,  the  Creatine+Spermine+Choline  to  Citrate  ratio, 
(Cr+Spm+Cho)/Cit,  can  be  measured  by  calculating  the  vol¬ 
ume  integrals  of  the  diagonal  and  cross  peaks  between 
2.9  and  3.3  ppm  for  Cr+Spm+Cho,  and  between  2.4  and 
2.8  ppm  for  Cit. 

RESULTS 

Phantom 

There  are  many  ways  to  view  the  4D  dataset.  One  simple 
way  is  to  view  each  2D  voxel  as  containing  an  individual 
2D  spectra.  Although  intuitive,  it  can  be  cumbersome  to 
view  all  spectra  at  once,  in  the  case  of  a  16  x  16  imaging 
grid  there  are  256  different  spectra.  For  simplicity,  a  cho¬ 
sen  peak  can  be  overlaid  on  top  of  a  localization  image  to 
show  the  spatial  distribution  of  that  peak,  as  shown  with 
the  strongly  coupled  (AB)  citrate  multiplet  (32,33)  (cen¬ 
tered  at  F2  =  2.6  ppm)  in  Fig.  2a.  Although  there  seems  to 
be  some  leakage  from  the  point  spread  function  as  well  as 
chemical  shift  artifacts,  the  peaks  with  the  brightest  inten¬ 
sity  appear  within  the  PRESS  excitation  volume  marked  by 
the  white  box. 

The  effect  of  NUS  leads  to  incoherent  aliasing  as  shown 
in  Fig.  2b  with  data  that  has  been  filled  simply  using 
zeros.  To  simplify  viewing  the  data,  one  column  (marked 
by  the  dotted  yellow  box  in  2a)  was  selected  to  display  the 
incoherent  aliasing  along  the  y-direction,  and  one  region 
centered  around  the  citrate  multiplet  (centered  at  F2  = 
2.6 ppm)  is  shown  to  display  the  incoherent  aliasing  in 
the  F1  dimension.  For  50%  of  the  data,  there  is  notice¬ 
able  leakage  of  the  citrate  multiplet  in  the  y-direction.  The 
peaks  from  the  fully  sampled  data  shown  in  2a  are  better 
resolved  than  the  under-sampled  peaks  in  Fig.  2b  (first  col¬ 
umn  in  each  pair)  and  still  maintain  the  same  shape  across 
the  PRESS  excitation  volume.  The  same  can  be  said  along 
the  F1  dimension  where  signal  leakage  is  taking  place  due 
to  the  NUS  in  Fig.  2b  (second  column  in  each  pair).  As 
the  number  of  sampled  points  is  reduced,  the  degree  of 
incoherent  aliasing  along  both  the  y  and  F1  dimensions 
increases.  At  16%,  most  of  the  original  shape  has  been 


heavily  degraded  as  the  signal  now  appears  to  be  incoher¬ 
ently  smeared  across  both  y  and  F1  dimensions.  By  10%, 
there  is  little  resemblance  to  the  original  dataset. 

The  CS  reconstruction  of  the  undersampled  phantom 
datasets  along  with  the  normalized  root  mean  square  error 
(RMSE),  calculated  with  respect  to  the  fully  sampled  data, 
is  shown  in  Fig.  3.  The  same  column  and  citrate  region 
in  Fig.  2  are  displayed  in  3a.  It  can  be  seen  visually  and 
from  the  low  RMSE  values  that  the  CS  reconstruction 
successfully  cleans  up  the  incoherent  aliasing  produced  by 
the  NUS  when  Figs.  2  and  3  are  compared.  For  as  little  as 
20%  of  the  original  sample  points,  the  spatial  reconstruc¬ 
tion  of  the  citrate  multiplet  is  very  close  to  the  original  fully 
sampled  data  (yellow  box  in  Fig.  2a).  At  16%  of  the  original 
data,  it  can  be  seen  that  while  the  spatial  distribution  along 
the  y-direction  is  still  consistent  with  the  fully  sampled 
data,  the  shape  of  the  peaks  starts  to  deviate  slightly.  The 
same  can  be  seen  for  the  reconstruction  of  the  F1  dimen¬ 
sion  where  the  algorithm  performs  well  in  removing  the 
incoherent  aliasing.  Again,  at  16%  of  the  original  sample 
points,  there  begins  to  be  slight  deviations  in  the  shape 
of  the  peak  (reflected  by  the  notably  higher  RMSE)  when 
compared  to  datasets  with  higher  percentages  of  the  origi¬ 
nal  sample  points.  At  10%  of  the  data,  it  can  be  seen  that 
while  a  lot  of  the  incoherent  aliasing  has  been  removed,  the 
final  reconstruction  has  noticeable  visual  deviations  from 
the  original  dataset,  along  with  an  RMSE  above  1%. 

Figure  3b  is  CS  reconstructed  spectra  from  only  20%  of 
the  original  data,  taken  from  the  same  location  as  the  exam¬ 
ple  spectra  in  Fig.  lb.  While  there  are  slight  differences  in 
noise  and  contour  levels,  the  same  basic  features  present  in 
the  fully  sampled  spectra  are  clearly  visible  in  the  recon¬ 
structed  spectra.  In  particular  is  the  faithful  reproduction 
of  the  citrate  multiplet  in  which  the  /-resolved  peaks  are 
nicely  resolved. 

In  Vivo 

Figure  4  shows  the  results  of  the  CS  reconstruction  of 
an  undersampled  in  vivo  scan  of  a  29-year-old  healthy 
prostate.  The  mask  that  was  used  to  sample  25%  of  the  kyt\ 
plane  is  shown  in  Fig.  4.  The  spatial  distribution  of  the  /- 
resolved  citrate  multiplet  at  2.6  ppm  is  shown  in  Fig.  4b.  As 
in  the  phantom,  the  peaks  are  localized  within  the  PRESS 
excitation  volume  (white  box)  with  some  slight  leakage 
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FIG.  3.  a:  Reconstructions  of  the  retrospectively  undersampled  datasets  shown  in  Fig.  2.  The  same  columns  in  Fig.  2  are  displayed  to  show 
the  effect  of  the  reconstruction  algorithms  along  with  the  normalized  RMSE  values,  b:  A  reconstructed  spectrum  from  the  same  voxel  in  Fig. 
1  b  with  only  20%  of  the  original  data  from  the  same  voxel  as  the  sample  spectrum  in  Fig.  1  b  and  contains  the  same  annotated  metabolites. 


due  to  the  point  spread  function  as  well  as  chemical  shift 
artifacts.  An  individual  /-resolved  spectrum  from  a  select 
2  cm3  voxel  is  shown  in  Fig.  4c.  The  extracted  spectra  have 
considerably  more  noise  than  the  phantom  scans,  but  the 
main  metabolites,  citrate,  creatine,  spermine,  choline,  etc. 
can  be  detected.  For  the  spectrum  shown  in  Fig.  4b,  the 
(Cr+Spm+Cho)/Cit  was  calculated  to  be  0.395.  The  mean 
(Cr+Spm+Cho)/Cit  value  for  all  six  volunteers  from  a  sim¬ 
ilar  location  was  measured  to  be  0.451  with  a  coefficient 
of  variance  of  19%,  and  is  consistent  with  what  has  been 
reported  in  the  literature  (34,35).  To  demonstrate  the  repro¬ 
ducibility  between  all  the  scans,  the  reconstructed  citrate 
multiplet  from  the  left  peripheral  zone  of  each  volunteer  is 
shown  in  Fig.  5. 

DISCUSSION 

The  application  of  CS  requires  that  the  fully  sampled 
dataset  be  sparse  in  some  transform  domain  and  that  the 
NUS  pattern  be  such  that  any  aliasing  be  incoherent.  Since 
the  undersampling  performed  in  the  kyt\  plane  of  the  4D 


EP-JRESI  data  was  shown  in  Fig.  lc  to  have  a  sparse  rep¬ 
resentation  using  finite  differences,  it  is  suitable  to  use  TV 
minimization  in  the  CS  reconstruction.  By  randomly  sam¬ 
pling  in  the  kyt\  plane  according  to  a  probability  density 
(Eq.  3),  the  data  yielded  incoherent  artifacts  that  spread 
across  each  NUS  dimension  as  seen  in  Fig.  2b.  The  inco¬ 
herent  artifacts  have  the  appearance  of  noisy  data  and 
get  progressively  worse  as  fewer  samples  are  used  in  the 
reconstruction.  While  the  incoherent  artifacts  are  not  tech¬ 
nically  “noise”  in  the  conventional  sense  in  that  they 
are  not  random,  they  may  be  removed  by  “de-noising” 
algorithms.  Total  variation  has  previously  been  used  as  a 
method  for  noise  removal  in  images  (36)  and  thus  serves 
as  a  suitable  objective  function  for  the  CS  reconstruction. 
The  reconstruction  can  be  seen  to  act  as  a  smoothing  algo¬ 
rithm  that  effectively  “de-noises”  the  data  so  long  as  it  is 
consistent  with  the  collected  data.  Such  data-consistent  de- 
noising  effectively  fills  in  the  missing  data-points  while 
maintaining  the  fidelity  of  the  originally  sampled  data. 

As  can  be  seen  in  Fig.  3,  the  CS  reconstruction  success¬ 
fully  cleans  up  the  incoherent  artifacts  for  undersampled 


FIG.  4.  Reconstruction  of  an  undersampled  4D  EP-JRESI  in  vivo  prostate  scan  with  only  25%  of  the  samples  as  required  by  the  Nyquist- 
Shannon  criterion  showing  (a)  the  mask  used  to  undersample  the  kyt i  plane  where  the  white  points  indicate  those  that  were  sampled,  (b)  the 
multivoxel  spatial  distribution  of  the  citrate  multiplet  at  2.6  ppm  overlaid  on  top  of  the  T2-weighted  axial  MRI  with  the  white  box  indicating 
the  PRESS  localization,  and  (c)  a  J-resolved  spectrum  from  a  select  voxel  in  the  reconstruction. 


1504 


Furuyama  et  al. 


F2  (ppm) 


2.7  2.65  2.6 

F2  (ppm) 


F2  (ppm) 


2.7  2.65  2.6" 

F2  (ppm) 


F2  (ppm) 


FIG.  5.  a-f:  Expanded  view  of  the  citrate  multiplet  from  a  select  voxel  in  the  left  peripheral  zone  from  each  of  the  six  healthy  male  volunteers. 


data.  For  the  undersampled  data  containing  more  than  20% 
of  the  fully  sampled  data,  the  reconstruction  produces  a 
high  fidelity  reproduction  of  the  original  data,  with  RMSE 
values  all  below  0.6%.  The  performance  of  the  reconstruc¬ 
tion  algorithm  begins  to  degrade  with  fewer  sample  points, 
in  this  case,  once  the  data  drops  below  20%.  Despite  the 
performance  drop,  it  is  still  remarkable  to  see  that  even  at 
10%  of  the  data,  the  CS  reconstruction  does  a  decent  job 
of  cleaning  up  the  incoherent  artifacts,  with  an  RMSE  of 
1.3%,  resulting  in  a  yF1  profile  that  is  similar  to  the  orig¬ 
inal  fully  sampled  data.  This  observation  can  be  useful  if 
a  high-fidelity  reconstruction  is  not  required,  and  the  data 
can  be  undersampled  to  a  higher  degree  to  save  time.  As 
expected,  the  RMSE  values  increase  with  the  acceleration 
factor  demonstrating  the  tradeoff  between  acceleration  and 
accuracy  of  the  reconstruction. 

The  undersampled  in  vivo  results  demonstrate  the  fea¬ 
sibility  in  reducing  scan  times  in  multidimensional  spec¬ 
troscopic  imaging  sequences  by  means  of  nonuniformly 
sampling  the  kyt 1  plane  and  CS  reconstruction.  While  col¬ 
lecting  only  25%  of  the  data  results  in  an  acceleration  factor 
of  4,  two  averages  were  acquired  for  SNR  purposes,  thus 
resulting  in  only  a  2-fold  acceleration.  Alternatively,  col¬ 
lecting  50%  of  the  data  with  only  1  average  would  result 
in  the  same  acceleration  factor  as  well  as  the  same  SNR. 
However,  without  sufficient  SNR,  many  of  the  lower  con¬ 
centration  metabolites  run  the  risk  of  being  de-noised,  or 
wiped  out  completely  by  the  reconstruction,  as  the  data 
fidelity  term  in  the  objective  function  allows  for  the  mod¬ 
ification  of  features  on  the  order  of  the  noise-level.  Since 
the  phantom  reconstruction  with  25%  of  the  data  was  com¬ 
parable  to  the  reconstruction  with  50%  of  the  data,  only 
25%  of  the  in  vivo  data  was  collected  with  two  averages  in 
an  attempt  to  boost  the  SNR  of  the  metabolites  with  lower 
sensitivity.  The  minimum  SNR  for  reconstruction,  which 


depends  on  the  overall  sparsity  of  a  particular  signal  as 
well  as  the  extent  of  the  incoherent  aliasing  is  an  ongoing 
field  of  research  and  will  be  the  subject  of  future  work  (37). 
Given  improved  sensitivity  (such  as  better  shimming,  more 
sensitive  coils,  and/or  higher  B0),  the  SNR  would  be  suffi¬ 
cient  without  having  to  average,  yielding  a  full  acceleration 
factor  of  4. 

As  the  4D  in  vivo  dataset  was  actually  undersampled 
during  the  scan,  there  is  no  direct  means  to  compare  with  a 
fully  sampled  dataset  as  was  the  case  with  the  reconstruc¬ 
tions  on  the  phantom  dataset.  The  reconstructed  in  vivo 
spectra  and  profiles  recorded  in  the  6  healthy  controls  show 
similar  quality  to  previously  reported  fully  sampled  EP- 
JRESI  spectra  acquired  in  the  prostate  of  healthy  volunteers 
(9).  The  citrate  multiplet  is  well  reconstructed  showing 
a  nicely  /-resolved  AB-type  structure.  The  other  major 
metabolites  are  visible  as  well,  in  which  their  calculated 
ratios  are  in  agreement  with  the  literature  and  consistent 
throughout  all  the  volunteer  scans. 

This  pilot  study  made  use  of  randomly  generated  sam¬ 
pling  masks  according  to  Eq.  3  to  selectively  delete  points 
in  the  fully  sampled  dataset  to  simulate  undersampling. 
The  same  type  of  randomly  generated  sampling  mask  was 
used  on  the  MRI  scanner  to  collect  an  undersampled  in 
vivo  dataset.  While  this  technique  satisfies  the  incoherent 
aliasing  criterion  in  CS,  it  is  not  necessarily  the  optimal 
sampling  scheme.  Not  all  sampled  points  have  equal  con¬ 
tribution  to  the  overall  structure  of  the  data,  and  it  can 
be  seen  that  the  collection  of  certain  points  can  be  more 
critical  than  other  points.  That  is  not  to  say  that  certain 
points  are  unimportant,  but  rather  the  collection  of  certain 
key  points  can  be  critical  such  as  the  point  at  ky  =  0  and 
fi  =  0.  Likewise,  it  is  possible  that  for  some  applications, 
one  of  the  dimensions  in  the  yF1  plane  may  be  more  sparse 
than  the  other,  and  so  the  values  for  a  and  b  (which  were 
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set  to  be  equal  in  this  study)  can  be  adjusted  accordingly. 
This  may  be  an  even  larger  concern  for  the  echo-planar 
correlated  spectroscopic  imaging  sequence  since  the  cross 
peak  signals  have  a  /-dependent  maximum  intensity  in 
the  t\  dimension  and  thus  cannot  be  modeled  as  simple 
decaying  exponentials  like  in  Eq.  3.  As  a  result,  deter¬ 
mining  the  optimal  sampling  masks  for  both  EP-JRESI  as 
well  as  echo-planar  correlated  spectroscopic  imaging  is 
the  subject  of  further  investigation  and  will  be  reported 
elsewhere. 

CONCLUSION 

In  the  acquisition  of  a  4D  spectroscopic  imaging  sequence 
(2D  spatial,  2D  spectral),  an  EPSI  readout  can  be  used  to 
simultaneously  accelerate  the  collection  of  one  spatial  and 
one  spectral  dimension.  This  is  the  first  in  vivo  study  to 
show  that  CS  can  be  used  to  simultaneously  accelerate  the 
acquisition  of  the  remaining  spatial  and  spectral  dimen¬ 
sions.  Two  techniques  have  thus  been  combined,  both  of 
which  simultaneously  accelerate  the  collection  of  one  spa¬ 
tial  and  one  spectral  dimension.  Such  acceleration  has  the 
potential  to  bring  4D  spectroscopic  imaging  (2D  spatial  and 
2D  spectral)  scan  times  well  under  the  coveted  10  min  bar¬ 
rier,  allowing  for  the  power  of  2D  in  vivo  MRS  to  become  a 
clinical  reality. 
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